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Silvia Polla, Philip Verhagen 


Introduction 


Within the framework of the Excellence Cluster Topoi, a fruitful interdisciplinary debate on 
space and movement over the long term has developed. The workshop “Computational ap- 
proaches to movement in archaeology” (organized on January 6 2011 in Berlin) tackled 
questions related to space and movement in the framework of computational archaeology, 
landscape archaeology, historical geography and archaeological theory. The current volume, 
which is the product of this meeting, brings together contributions that show how the study 
of settlement patterns and movement has been dramatically transformed by the use of spa- 
tial technology (GIS), in particular Cost Distance and Least Cost Paths (LCP) Analysis. The 
term “least cost path" is somewhat deceptive, however, since it is not just the costs of move- 
ment, but also the benefits of moving to a particular location that influenced the routes 
chosen and created. 

Archaeological theories about the way people moved in the landscape, and how they 
created and maintained paths and communication networks are often based on relatively 
abstract notions. For example, several papers in the current volume indicate that visibil- 
ity may have been an important factor (co-)determining movement and path creation in 
the landscape. However, the exact parameters involved, and how they influenced the 
routes chosen, are largely within the realm of speculation. Computer-based modelling 
can be seen as a sophisticated approach to speculation. It allows us to experiment with 
the possible parameters involved, change the values and weights of each and inspect 
the outcome to see whether it conforms to our initial expectations and if it in some way 
fits the actual archaeological evidence. Most importantly, computer-based models are 
explicit: since all assumptions are laid out in detail, we can study the consequences of 
changing them, and the models can be replicated. It is through modelling that different 
scenarios can be explored and compared to real-world outcomes. Computer-based mod- 
els are therefore in essence heuristic tools that can help to develop theory and interpre- 
tation. 

The papers in this volume all share this perspective: path modelling is used not so 
much as a tool to reconstruct ancient routes, but rather to help the interpretation of existing 
or hypothesized pathways or to develop new theories about prehistoric movement. They 
discuss the theory and methodology of path modelling and use case studies to show the po- 
tential of their approaches for understanding factors and effects in landscape formation 
and transformation in the long term. The issues involved are complex, however, and the 
theoretical and methodological debate focuses on four main issues: 

I the extraction of indicators of accessibility and potential movement at different scales; 
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2. model building procedures involving the choice of parameters and the definition of 
the cost functions; 
3. theevaluation of the created models (best fit, deviations, required accuracy); 

and network analysis approaches, especially concerning specific cultural contexts like 

Roman infrastructures and coherent road systems. 

The interpretation ofthe generated models of networks and their relation to settlement pat- 
terns is concerned with time, scale and mobility issues, such as the type and purposes of 
movement and infrastructures, attractors of and barriers to movement. Continuity and 
change in settlement patterns can be studied by addressing the chronology and stability of 
paths and communication networks through a landscape archaeology approach. 

The contribution by Dimitrij Mlekuz discusses the applicability of the concepts of time 
geography to issues of past movement. The theoretical framework he sketches is greatly in- 
fluenced by postmodernist philosophy and by the notion of time in landscape archaeology 
studies. Movement is fluid, and constantly re-creates the material world. The dynamic na- 
ture of movement is problematic from a technical point of view, however. GIS-based map- 
ping only offers a detached view, presenting fixed snapshots in time. Mlekuz argues that we 
can escape from this position by regarding computer-based models as practices that allow 
us to explore the connections between static representations and the flows and movements 
of people and objects in the past. Time geography is then used as a conceptual structure for 
creating "potential path areas," basically the area that can be reached within a certain travel 
time. This method creates what he calls a topography of accessibility. While not completely 
new from a technical point of view, the theoretical framework and the archaeological possi- 
bilities of interpretation that this approach offers are promising, since it relates standard 
cost surface calculations to issues of scale of movement (long distance versus daily practice) 
and strategic considerations in the placement of settlements. 

The paper by Gary Lock, Mariza Kormann and John Pouncett focuses on the possibil- 
ities of integrating and better understanding the role of visibility in creating prehistoric 
pathways. Their paper is very much inspired by a post-processual, phenomenological the- 
oretical framework, and is noteworthy for the attempt that is made to understand in detail 
what it actually means to include visibility in least cost path modelling, and how quantitative 
and qualitative approaches can be integrated in both method and theory. Lock et al. have 
developed a sophisticated multi-scalar approach to path modelling, based on viewsheds that 
take into account the distance to the horizon. Futhermore, the weight attributed to visibility 
as a factor influencing movement is modified according to the possible goal of movement: 
do people want to see around them, or do they want to avoid being seen? The resulting paths 
and interpretations, using a case study of Neolithic long barrows in England, show that path 
models should primarily be seen as tools for exploring the connections between the grave 
monuments and movement at what could be called a *phenomenological scale." 

Lock et al.'s paper is strongly connected in methodology to the contribution of Élise 
Fovet and Klemen Zaksek. They use a different approach to modelling visibility, using the 
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average angle of view rather than the distance to the horizon, but the aim is very comparable: 
to find out how visibility influences the location of the modelled paths. In their case, the 
interest is in reconstructing and interpreting the Roman road network in the Languedoc 
(France). In contrast to the first two papers, however, they are not particularly interested 
in the phenomenological aspects of movement, but much more in the long-term evolution 
of settlement patterns and communication networks, and the relations between settlement 
location and routes. 

The paper by Philip Verhagen, Silvia Polla and lan Frommer focuses more on metho- 
dological aspects of route modelling. In a similar vein to the previous two chapters, they 
try to reconstruct a Roman road section in Cappadocia (central Turkey) using visibility as a 
potential factor influencing road construction. To judge from the modelling results, how- 
ever, this does not seem to have had a major influence. The Roman road network in the area 
was completely restructured during the Byzantine period. Instead of long-distance connec- 
tions, shorter distance roads were constructed. One of the features of this more fine- 
grained network of connections is that junctions are found in places where no settlement is 
present. Modelling the location of these junctions is not possible with standard GIS tools, 
however. Instead, an algorithm for creating junctions (so-called Steiner points) was experi- 
mented with, and resulted in modelled paths that are closer to the supposed location of the 
Byzantine roads than standard least cost path calculations. 

The last chapter, by Patricia Murrieta Flores, examines in detail the role of navigation 
in movement, and specifically in the patterns of transhumance in the Copper Age of Anda- 
lusia (Spain). The presence of megalithic funerary monuments in the mountainous area of 
the Sierra Morena is taken as the point of departure for examining the possible role ofthese 
monuments as territorial markers that structured the seasonal movement of people and 
animal herds. The statistical analysis of the visibility of these monuments, in conjunction 
with the modelled natural corridors of movement, indicates that the megaliths may not 
have been navigational markers per se, but certainly marked vistas and may have been con- 
nected to a symbolic meaning as places of the ancestors, visited annually in the transhu- 
mance cycle. 

As is evident from the contributions in this volume, path modelling is a subject that 
has proved to be one of the more challenging research topics within spatial analysis in ar- 
chaeology. Considerable progress has been made over the past decade towards a better 
understanding of the mechanics of movement and in improving the methods to model 
this, and the archaeological implications of using path modelling are the subject of lively 
debate. So far the actual impact of path modelling in (landscape) archaeological research is 
still relatively limited, and its use has probably raised more questions than it has provided 
answers. In itself, this is not a development that we should be overly worried about. The use 
of new techniques and novel theoretical concepts will always open up new horizons as well 
as blind alleys and, from the successes and failures, fresh perspectives on the nature of 
(pre-) historic movement patterns, path finding and the development of communication 
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networks will result. However, as the various contributions in this volume show, path mod- 

elling needs to be combined with well-framed archaeological research questions and sup- 

ported by additional analyses to be used to maximum effect. And as long as path modelling 

still demands a relatively large investment in learning to master the technique and deal 

with the methodological problems, it will not become standard fare for spatial analysis — 

a Characteristic it shares with the related fields of visibility analysis and network analysis. 

A potential research agenda for the next five to ten years should therefore focus on the fol- 

lowing major issues: 

- improving the efficiency, flexibility and general ease of use of the available software 
solutions; 

— developing methods and tools for better integration of path modelling with related 
fields like visibility analysis and network analysis; 

— research into the mechanics of non-pedestrian movement, and the effects of other fac- 
tors than slope on movement speed; 

— developing better methods for testing the results of path modelling using archaeologi- 
cal data and theoretical concepts; 

- and, finally, translating archaeological theory concerning movement and communi- 
cation into terms that will allow for easier modelling. 

We hope that the current volume may serve as a source of inspiration for archaeologists to 

follow up on these issues. We are convinced that path modelling will continue to be an ex- 

citing and fruitful area of archaeological research for years to come. 
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Dimitrij Mlekuz* 


Exploring the topography of movement 


1 Introduction 


Michel de Certeau, at the beginning of his groundbreaking essay “Walking in the City,” 
climbs to the top of World Trade Centre to enjoy the view of New York laid out below 
(de Certeau 1984). However, he notes that “to be lifted to the summit of the World Trade 
Centre is to be lifted out of the city's grasp." This is a privileged view, but one completely 
detached from the city below. 

A city can exist only if there are people in it. The ordinary practitioners of the city are 
the pedestrians who make use of the spaces to walk, and they bring life to the city. Walking 
is “an elementary form of this experience of the city” (de Certeau 1984, 94) and integral to 
this experience are walkers, “wandersmänner, whose bodies follow the thicks and thins of an 
urban ‘text’ they write without being able to read it" (de Certeau 1984, 93). 

Walking in the city turns out to have its own logic, or, as de Certeau puts it, its own rhet- 
oric. The walker individuates and makes ambiguous the "legible" order given to cities by 
planners, which can be seen on maps or from above. 

The same can be said for landscapes. Movement constitutes landscape; it is an essential 
act of weaving places into the web of landscape, making it real for the people themselves. 

People also carry objects, moving them into new positional and relational contexts 
with other things, and they create new material encounters and new material traces 
through movement. 

All of these parts are moving: movement suggests a continual process of becoming. 
Places are not so much defined by their location, their boundaries or shape, but through the 
flows and convergences that flow through them. Places are gatherings that are continuously 
on the move. 

Movement is thus continuously generative and, if we want to understand it archae- 
ologically, we need to develop approaches that can cope with the fluid conditions of these 
generative practices. 

Representing movement, pinning it down on maps, has the effect of arresting move- 
ment outside the flows of its temporal and spatial contexts. Yet our discipline requires us to 
give movement an image. The common GIS approaches to movement are problematic, as 
they show not the emergent aspects of movement, but a detached view, de Certeau's view 
from above. 


* University of Ljubljana. 
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But that does not mean that our attempts to tackle movement in the past using com- 
putational approaches are missing the point. Attempts to represent movement should not 
be seen as points that record and communicate information about fluid and generative 
practices. Instead they should be seen as activities that contribute to the constitution of these 
fluid practices. Attention should shift from the image and its fidelity or treachery, towards 
practical relations and processes, which actively produce truths (Wickstead 2009, 254). 

Computational tools help to articulate and express what we mean, because they are 
tools that help us to intertwine our practices of modelling with past movements. As Aldred 
and Sekedat (2011a) write, “an archaeology of movement should not exclusively focus on 
the captured image of representation, but also contend with the interdependent and emerg- 
ent meanings that arise during intervention." 

This paper is thus an attempt to approach movement from a computational perspec- 
tive, creating tools that can help us cope with the fluid and emergent nature of movement, 
but instead of producing fixed representation in maps, we should see them as practices — 
open-ended, processual, explorative activities - which do much more than resemble: they 
take the places of the original situation and thus enable us to intertwine our practice of 
modelling with the flows and movements of the past. 


2 Time geography 


In movement, issues oftemporality, spatiality and practices are dialectically woven together 
and inseparable. All these issues are addressed holistically by time geography, as developed 
by Torsten Hägerstrand and his associates at the University of Lund in the 1970s (Háger- 
strand 1970; 1973; 1975; 1976; Pred 1977). Time geography is a powerful conceptual frame- 
work for understanding human practices in both time and space, and, in particular, con- 
straints and trade-offs in the allocation of limited time among practices in space. Time 
geography focuses on constraints imposed by the corporeality of human existence in time 
and space, such as the fact that no-one can be in two places at the same time, that everyone 
has only limited ability to undertake more than one task at a time, that movement is time 
consuming and “the fact that every situation is inevitably rooted in past situations" (Háger- 
strand 1975; Pred 1977). 

Time geography has developed a series of analytical instruments to investigate how 
certain types of constraints have a structuring effect on everyday life. The most fundamen- 
tal concept of time-geography is time-space, or the time-space aquarium, a bounded three- 
dimensional space where width and length describe two-dimensional space, and depth 
stands for time. The time-depth of the aquarium may vary, from daily or yearly, or a bio- 
graphical, lifelong time scale. 

Any physical body (a thing, substance, animal or human person) within the time-space 
aquarium can be described by its time-space path or trajectory, which describes its, his or 
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her position and movement within the time and space. Thus, in time geography, time and 
space are inseparable. 

However, much more interesting than isolated paths of agents are the time-space pat- 
terns, that emerge from interactions and relations between agents. This “choreography of 
existence” (Pred 1977) emerges when the movement of agents on their daily practices 
forges a web of relations between them that serve as further constraints. Thus, for example, 
to join in a shared activity, two or more agents form a group and their time-space paths join 
at some place. The availability of material, tools and equipment is a further constraint on 
where this place could be. Through such an activity, systems emerge: these are stable, 
multi-scalar spatiotemporal patterns formed from the intertwined allocation of time among 
practices in space. Thus, in time geography, movement is not understood in isolation but 
instead as a component in a wider pattern of interaction, as “the glue in connecting” as 
Aldred and Sekedat (2011b) put it. 

The landscape - the relation between things where agents move around actively — 
interacts with the capabilities of agents and thus plays an active role in the shaping of the 
paths. Hägerstrand (1976, 332) calls it the “principle of togetherness;" this is not, he writes, 
^just resting together. It is also movement and encounter ... because what is all the time 
resting, moving and encountering is ... humans, plants, animals and things all at once." 
Thus movement is the intertwining of many lines that run concurrently, making it a gen- 
erative and recursive practice. 


2. Time geography and archaeology 


Matters become more complicated if we try to apply time geography to archaeological prob- 
lems. Due to the fragmented nature and resolution of the archaeological record, the direct 
applicability of these concepts to archaeology is extremely limited. Except in outstanding 
cases, it is impossible to trace a time-space path of persons in the past. However, Háger- 
strand and other time-geographers point to the regularities in movement when people 
carry out their everyday tasks or projects, which emerge as an activity system. Activity sys- 
tems are patterns of daily, seasonal or lifetime routine movements, interactions and activ- 
ities by agents. Practices can thus be understood as material bodies of work that have 
gained enough stability over time, through the establishment of corporeal routines and 
specialized devices, to reproduce themselves (Thrift 2007, 8). 

Different persons tend, for various reasons and in various temporal scales, to move 
regularly in similar paths in time and space, for example going to work in a field, travelling 
to market, herding sheep, etc. Social practices or routines are normally conducted at defi- 
nite places and at given moments (“stations,” in time geography speak) although some may 
extend over long periods of time. Thus social individuals tend to follow recurring time- 
space paths in their day-to-day lives. Routine practices are shaped by various forms of con- 
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straint, for example the need for regular sleep and food, the seasonal requirements of plants 
or animals, the availability of tools, the availability of other agents who form activity 
bundles, constraints imposed by the landscape itself, etc. 

The routine practices of course also change the places where they occur; they leave 
material traces that create new constraints in the process. In this way they have structuring 
effects on future practices that will occur at the station or in the vicinity. Every movement 
changes the relation between things (Fahlander 2003, 122-132). As Julian Thomas (1996, 
90) puts it: “[w]hile people often move in cyclical patterns in the course of routine activities, 
returning to the same location again and again ... the places ... are themselves continuously 
being physically altered and decaying, as well as continuously being re-evaluated and re-in- 
terpreted.” 

The archaeological record is of course not a result of acts derived from a single person, 
but the cumulative effect of practices, brought into new relations with our own practices of 
archaeological work. Agents on their daily routines have shaped and changed the land- 
scapes over a long time. The settlement then becomes the material residue of a bundle 
of special activities, which in turn has enabled other activities to occur in its vicinity, includ- 
ing the work of archeologists. Thus place is not an accretion of traces over time, but a com- 
plex set of spatial relations in a temporal setting, derived from the movement through it 
(Aldred / Sekedat 2011b). 

We can thus see the landscape as a network of places, which are constantly reproduced 
as “collapsed acts” of persons moving and performing routines in time and space. This is 
close to Tim Ingold's (1993) relation between the “taskscape” and landscape: “just as the 
landscape is an array of related features, so — by analogy - the taskscape is an array of related 
activities." 

This forces us to explicitly address the notion of time, which in this case is not just the 
simple linear succession of events, but is itself structured and shaped by the performance 
of tasks and produces different “temporalities,” which are “collapsed” into landscape. Thus 
“[IJandscape is time materialized. Or, better, landscape is time materializing: landscapes, 
like time, never stand still" (Bender 2002, 103). 

Landscapes are in a constant process of becoming (another temporality) through the 
daily, seasonal or biographical time-space “ballet” of actors. The routine practices in the 
landscape thus have structuring effects, which can be seen in the material record. As people 
move around in their daily practices, they make landscapes and make time go on. 
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3 Topography of accessibility 
3.1 Time-space prisms 


How can these fluid, generative, continuous, aspects of movement be tackled using GIS? In 
archaeological analyses, past movement is usually modelled using cost surface generation 
and the computation of “least cost" pathways between given locations in the landscape. 
A cost surface is a continuous surface in which each part of the landscape (i.e., raster cell) 
represents the cost, or “price” to reach that point from a source point. Cost is computed 
from the distance and "friction" or difficulty of moving across the raster cell is expressed 
either as time or caloric expenditure. Thus price can be expressed as an accumulated fric- 
tion (Bell / Lock 2000, 86). Least cost pathways between the source of cost surface and tar- 
get are traced by using some kind of optimizing algorithm (usually a derivate of Dijkstra's 
shortest path [1959] algorithm). The shape of the computed path is extremely sensitive to 
the algorithms used and so least cost pathways can diverge significantly. 

The typical applications of cost surface analyses in archaeology, for example analyses 
of site catchments, least cost paths etc., usually implicitly or explicitly invoke theoretical as- 
sumptions such as rational behavior and the least cost principle, as the name “least cost 
analysis" implies. 

Thus an “optimal path" is represented as a static map, line or series of points which 
connects the departure and destination points. It looks very precise, but its accuracy is ques- 
tionable. 

A more general use of *optimal paths," which is in many ways similar to approach out- 
lined in this paper, was attempted by Marcos Llobera (2000). His approach focuses on a 
general pattern of movement in the landscape rather than specific pathways, and he ex- 
plicitly approaches movement as a vital element of social landscapes. 

We argue that cost surface analyses can have a much more general analytical use in ar- 
chaeology. We turn for inspiration to time geography. 

The future paths of a past agent located somewhere in space can potentially branch off 
in a number of directions. However, they are limited by constraints working on them, such 
as maximal speed, availability of transport, etc. The time-space prism circumscribes a time- 
space volume where all trajectories must or can progress. Its shape is defined by constraints 
that operate on the individual. The prism thus discloses possibilities, and focuses on these 
rather than on behaviour per se. This renders the concept neutral and universally applicable, 
as it articulates a fundamental coherence in our world (Lentrop 2004). 

The time-space prism can exist between any pair of temporally adjacent points. In this 
case, there is an open temporal interval — “time budget,” in time geography parlance — dur- 
ing which the person can conduct travel and participate in activities. More generally, one 
or both of the origin or destination points may be unknown, or they may be different 
locations. Projecting the time-space prism onto the two-dimensional geographic plane de- 
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limits the potential path area, that is, the set of geographic locations that the person can oc- 
cupy during the available time budget. The potential path area is a binary map, where 
values indicate presence inside the projected time-space prism. 

The cost surface, where the cost represents the travelled time from the starting point, is 
an equivalent of the surface of the time-space prism with one known and one unknown point. 
We might call it the future surface, as it represents the time-space surface of maximal possible 
extent of an agent at a given time after she departs from the starting point. In the case that the 
point represents the known end-point, the cost surface represents the past surface, or possible 
locations at a given time before she reached the destination point (Miller 2005). 

Most of the widely used cost surface analyses in archaeology, for example the delimi- 
tation of isochronous surfaces around sites (Gaffney and Stančič 1993), are therefore 
potential path areas of prisms, where agents return to their destination points within a 
given time budget. These analyses can of course be produced using ready-made capabilities 
offered by standard GIS packages. 

Potential path areas are the most general way of describing movement. They focus 
only on constraints that operate upon the moving actor. In this way they show where 


DIMITRIJ MLEKUZ 


movement could occur, given the constraints. It is extremely imprecise, but also very ac- 
curate. 

However, there are also some more exciting uses of potential path areas, which can 
yield a more complex or deeper understanding of the social processes in the landscape and 
tell us more about how people move around and perceive the landscape. 

We can combine potential path areas in a new map by adding different potential path 
areas from a set of starting points. In this way we get cumulative path areas. The value of 
every cell in the cumulative path area map indicates from how many starting points it is ac- 
cessible within the allocated time budget. By modifying the time budget we can observe the 
effect of it on the structure of possible interaction in the landscape. Thus we can measure 
the accessibility of parts of the landscape from the different sites or bundles of activities. Its 
shows how the places change the structure of the landscape, and can identify the places or 
areas where interaction is facilitated by the landscape, where bundles of activities can occur 
and where not, or where a larger time budget is needed to form a bundle. It can also help 
the identification of *non-events," or better, the barriers that prevent certain types of activity 
from occurring (Pred 1977, 210). 


3.2 Potential path field 


The idea of cumulative path area can be expanded to include every point in the landscape. 
The potential path field is thus a sum of potential path areas within a given time budget from 
every location in the landscape. Thus the potential path field is computed by adding poten- 
tial path areas from every location in the landscape within the given time budget. The result 
is a continuous raster surface, where the value of each cell informs us from how many 
points within a given time budget it is accessible (Fig. 1). Every point will have a value of at 
least 1, as it is accessible from itself, but usually a much higher value. Points with higher 
values are accessible from more locations within the same time budgets, or to turn it 
around, they can access more locations than points with a lower count. 

Potential path fields are computationally extremely intensive. The results also suffer 
badly from the boundary effect, which means that the values at the edges of study area do 
not provide real values. 

Potential path fields are at first difficult to interpret. They give a cumulative overview of 
accessibility, and do not describe specific acts of movement, or movement from a specific 
location. They therefore yield a new idea of landscape, based on its accessibility. We can also 
call it the topography of accessibility. It thus identifies areas that are more accessible and so 
shows the structuring potential of the landscape itself on the movement and practices in 
landscape. 

By varying the time budget we can explore the relation between time and the structur- 
ing effect of the landscape (Fig. 2). The topography of accessibility is highly dependent on 
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time budget. With a limited time budget, absolute values are generally small, as one can ac- 
cess only a limited number of locations from each spot. Local differences thus have a huge 
effect on the topography of accessibility. Differences between close locations are therefore 
sharply pronounced. Topography of accessibility with a small time budget shows fine- 
grained differences in accessibility. As time budget increases, the topography of accessibil- 
ity becomes absolutely higher, but more and more smoothed out and less grainy, as close lo- 
cations become more and more similar in terms of accessibility. We can observe the emerg- 
ence of large-scale patterns, such as accessibility basins, barriers and corridors. 

Varying the available time budget allows us to explore how time spent performing 
practices changes the way people interact and access the landscape. Thus locations with 
good accessibility on a small time budget do not necessarily have good accessibility when 
the time budget is increased. The landscape activities that require more time (for example, 
going on pilgrimage, transhumance, trade ...) structure or define landscape in completely 
different ways from practices of short duration (for example, daily activities, going to the 
fields, herding sheep, carrying substances to a house ...). Landscape, at least from the per- 
spective of accessibility, is therefore constructed from practices that are performed in it. 

Thus, in addition to busier areas in the landscape, we should also expect to find non- 
places and back-areas. We can measure how accessible parts of the landscape are, where in- 
teraction is facilitated by the landscape, where bundles of activities can occur and where 
not, or where a larger time budget is needed to form a bundle. In this way we can see land- 
scape as a vehicle for interaction, and the ways daily practices structure the long-term ma- 
terial record in the landscape. 

In short, the potential path field or topography of accessibility shows how the topogra- 
phy affords interaction and access on different time scales, or how landscape facilitates and 
is structured by practices. In this perspective, landscape becomes a very dynamic and fluid 
thing. Practices that differ in their term of duration or available time budget create different 
landscapes. 

The potential path field is in many ways similar to Marcos Llobera's (2000) topo- 
graphic accessibility index (or total path cost). His index is computed by summing all opti- 
mal paths from each cell within a given neigbourhood to a destination cell, for each location 
(cell) in the landscape. However, the potential path field, being developed from the perspec- 
tive of time-geography, explicitly references time as the main constraint that limits the 
movement and interaction in the landscape. Here, movement is approached in the most 
general way, through time-space prisms expressed as cost surfaces, making the potential 
path field a more general and robust index than Llobera's approach, which focuses on the 
specific paths. 
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Ocra pass 


Fig. 2 | Digital elevation model of the study area and topographies of accessibility for different time budgets. 
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3.3 Working with the topography of accessibility 


The potential path field or topography of accessibility measures the connectedness of each 
cell to other cells reachable within the time budget. Higher on the topography of accessibil- 
ity, more places become accessible within the given time budget: climbing the high ground 
one gets better access to the places within the time budget; staying on the high ground 
maintains good access to other places within the time budget, meaning that travel on the 
high ground is fast. Cells positioned high on the topography of accessibility are thus candi- 
dates for paths, corridors and routes. 

If we see the potential path field as the topography of accessibility, we can analyse the 
topography of accessibility using morphometric parameters such as slope, aspect and other 
features (Hengl / Reuters 2009). 

Thus the slope of topography of accessibility (or first derivate of the potential path 
field) gives us the change in speed. In steep areas travel slows down or speeds up, depend- 
ing on whether one is travelling up or down on the topography of accessibility. The steeper 
the slope, the larger is the change in speed. “Flat” areas indicate areas where speed of move- 
ment is constant (Figure 3). 

The aspect of the topography of accessibility is the direction where the change in ac- 
cessibility is largest. 

Using morphometric tools, the topography of accessibility can be subdivided into fea- 
tures such as pits, peaks, channels, ridges, passes and planes. The names of features suggest 
a geomorphological interpretation, but they must be understood in terms of rates of 
change of three orthogonal components of the potential path field. Pits are the locally least 
connected cells, while peaks are the locally best connected ones. Ridges connect points 
with good accessibility, making them natural corridors for movement. On the other hand, 
valleys are poorly connected, backwater areas. Saddle points are points between two ridges, 
points or areas where backwater areas can be crossed. In this way they can be seen as 
bottlenecks that funnel movement. Plains are areas where speed of movement is constant 
(Fig. 3). 

A prominence filter (Llobera 2001) reveals the difference of accessibility between the 
point and other points in the landscape within the arbitrary neighbourhood. The accessi- 
bility prominence is defined as the percentage of locations that are less accessible than 
the individual's location within a certain radius. In this way a prominence filter isolates 
parts of the landscape within the neighbourhood that are more accessible within the spe- 
cific radius, that is, locations that are better connected than their neigbours. Landscape is 
then structured into a network of well-connected locations and islands of backwater lo- 
cations (Fig. 3). 

All these analyses make it possible to explore how landscape is structured into 
pockets of backwaters and corridors of movement on different temporal and spatial 
scales. All these more or less distinct features, places or locales may be important for 
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Morphometric 
features 


Fig. 3 | Morphometric analyses of the topography of accessibility. 


understanding how people really move around the landscape and how they perform 
activities there. In a time geography perspective, such places and features can function 
as corridors, which channel the flow of movement, structure daily practices and enable 
bundles of activity to occur. 


4 Instead of a case study: Spatial strategies and tactics around the Ocra pass 
in prehistory 


How can these tools be used to understand past movement? Instead of a full-fledged case 


study we present a preliminary sketch that highlights the contributions of the new tools to 
the understanding of various practices in the landscape. We demonstrate how these tools 
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Fig. 4 | Morphometric features of the topography of accessibility of the study area on a 60 minute and a 5 minute 
time budget. 


can be used in a more nuanced and open-ended exploration of past practices in a land- 
scape.' 

Razdrto, or Ocra, as it was called in classical sources, is a pass around 597 m a.s.l. 
crossing the Dinaric Alps (Alpes Iuliae in classical sources) in western Slovenia. Ocra is the 
gateway between the Mediterranean and Central Europe. Here, two old long-distance chan- 
nels funnel: the Amber Route connecting Adriatic Sea with the Baltic Sea, and the route to 
the Black Sea, from where Argonauts came on their mythical journey from Colchis. 

In the Bronze and Iron Ages a series of hillforts was constructed in the region, some of 
them directly overlooking the pass. Roman transport across Ocra increased by the mid-2nd 
century BC, after the establishment of the colonia Aquileia in the Friuli Plain. In the Au- 
gustan period, a road was built, which connected Aquileia with the Roman town Emona 
(Horvat / Bavdek 2009). 

How can we approach the area of the Ocra pass in terms of mobility? Surely there is no 
need to tackle obvious paths and corridors. The pass, which *naturally" channels the move- 
ment, is very obvious, as are the very easily observed general corridors to and from the pass 
(Fig. 2, 3 and 4). 

The 60-minute potential path field clearly includes areas with long-distance, or long- 
term accessibility, areas which are very well connected with distant locations (Fig. 4). It 
shows the main corridors where people and things can flow relatively fast and unhindered 
by topography. This is accessibility in terms of fast movement across the landscape and not 
in the landscape. 


I For this study, a DEM (digital elevation model) ofthe 60.5 km? study area with the spatial resolution of 12.5 m 
was used. Friction surface was computed using Tobler's hiking formula (Tobler 1993) and expressed as minutes 
needed to cross the raster cell. All calculations were performed using scripts developed for GRASS GIS. 
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The establishment of hillforts around the main corridors is thus an element of a plan 
to control and channel the flow of people and material across the pass. The communities 
who built them created new spatial arrangements, which were made to control and domi- 
nate the passage of the material and people across the pass. 

This obvious, visible spatial arrangement can be identified with what de Certeau calls 
“strategy.” For de Certeau (1984), power is about territory and boundaries, in which the 
weapons of the strong are classification, delineation, division and control, while the weak 
use furtive movement, short cuts and routes (his term is “tactics”) to contest this spatial 
domination. 

Strategies are only available to subjects of “will and power,” so defined because of their 
access to a spatial or institutional location that allows them to objectify the rest of the social 
environment. Strategies thus invoke and actualize a schematic and stratified ordering of 
social reality. 

The position of hillforts around the obvious channels and corridors is thus a strategy, 
by which established centres of power are made concrete in hillforts in relation to the flows 
of material across the pass. 

As we decrease the time budget, to 5 minutes for example, a more surprising picture 
emerges. The landscape composed of major corridors fragments into a mosaic of small- 
scale paths and routes outside the grand corridors (Fig. 4). There is no longer a single, well- 
ordered landscape structured around the main corridors of movement, but a chaotic, unor- 
dered landscape. In De Certeau's terms, this is the space of “the tactic," an action that he 
defines as insinuating itself within the space of the other, worming its way into the territory 
of that which it seeks to subvert. It claims no space for itself, relying rather on time - “it is 
(de Certeau 1984, 


I» 


always on the watch for opportunities that must be seized “on the wing 

Thus we see a web of possible movements, crosscutting the structure ordered by inter- 
community politics and power, materialized in the hillforts. It is the place of daily action, of 
the practices of people doing their daily chores outside the domain of politics and large 
structures of power. 

These places possibly never materialized into structures like paths, as the gains from 
their small-scale victories are always discarded, *whatever it wins, it does not keep" (de Cer- 
teau 1984, xix). A tactic is to infiltrate into the other's place, always on the watch for oppor- 
tunities in time. For de Certeau, this is the central feature of everyday life, which *poaches" 
in countless ways from the property of others. People appropriate and reappropriate the 
strategic structures, in order to manipulate and enjoy them. The hillforts make sense only 
in relation to the large channels of movement, that is, movement in terms of long-distance 
communication, transport and exchange. When we get down to the tactics of daily prac- 
tices, however, the whole landscape, ordered and constructed by the logic of power and con- 
trol, fragments into the landscape of daily practices. The daily practice of walking subverts 
the logic of control and domination. 
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In this way we can explore the ways people resist the spatial discipline imposed by 
communities. This “imprint of walking” breaks down the logic of domination through the 
spatializing practices of everyday living. 

It is ultimately the people, on their daily tasks, who create landscape. As with de Cer- 
teau's walkers in the city, landscape is created by daily tasks in the landscape and not by 
large-scale plans and ideas. Landscape is not a plan laid down by strategies of control and 
domination, which forces people into predefined channels of movement. Ultimately it is 
people who trace or do not trace these corridors, using their tactics to avoid them or subvert 
them for their own purposes. 


5 Conclusions 


As people move through the landscape, so they define it: there is no landscape without 
people moving around on their tasks. People create the space to move through as they move 
through it, and their movement creates time as well. The real challenge facing archaeologi- 
cal GIS research lies in the development of concepts, tools and approaches that can enable 
us to analyse and understand the complexities of past human engagements with the world. 
Time geography offers a conceptual and methodological framework for conducting accessi- 
bility constraint analyses that can be implemented within GIS. 

Time geography focuses on time, space and practices, and places the emphasis on the 
significance of the material world and its constraints. Time geography is not limited to the 
mere description of courses of events in time and space, but emphasizes the importance of 
gaining a deeper understanding of fundamental social conditions and processes (Lenntrop 
2004, 223). The approach outlines necessary but not sufficient conditions for human inter- 
action. Thus it provides an approach rather than a full theory and requires the utilization of 
other knowledge and theoretical frameworks. We see this as its this major strength, as it 
does not come loaded with theoretical baggage, but can be employed as an analytical tool 
within many different theoretical approaches towards past landscapes. 

By focusing on the human interactions afforded by time and space, it allows us to ap- 
proach the landscape as an agent of interaction between people. A time-space framework 
can help us to understand the limits time and space impose on interaction between people. 
It tells us about possibilities and nothing more. It does not offer spatially precise but non- 
accurate results, such as are provided by least cost paths. 

Finally we want to stress that our tools do not produce non-problematic represen- 
tations of movement, as end-points that record and communicate information about move- 
ment: paths, roads, corridors ... We believe that those representations would ultimately be 
misleading and wrong. Instead, these tools should be seen as explorative activities, that do 
much more than just resemble: they take the places of the original situation and so enable 
us to intertwine our practice of modelling with the flows and movements of the past. The 
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difficult task of figuring out what happened within those possibilities is thus left to the 
interpretative skills of archaeologists. 
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Gary Lock, Mariza Kormann, John Pouncett* 


Visibility and movement: towards a GIS-based integrated approach 


1 Introduction: why visibility and movement? 


A continuing theme since the initial adoption of GIS technology by archaeologists has been 
the modelling of visibility and movement as ways of understanding past landscapes and the 
people who inhabited them. The underlying issues have been fully explored elsewhere (Lock 
/ Pouncett forthcoming) and here we briefly introduce only the main points relevant to the 
work presented in this paper. Modelling visibility and movement is usually justified through 
an interest in “perceiving the landscape" and an attempt to move GIS-based analysis away 
from "detached objectivism" (describing landscape as if viewing from a distance) towards 
"situated subjectivity" (trying to understand landscape from experiencing it directly). 

From the earliest publications on archaeological GIS, the potential of modelling vis- 
ibility and movement has been recognized. In his introduction to Interpreting Space, for 
example, Green talks of *perceiving landscape" and how *GIS provides such a means" 
(Green 1990, 6), while within the same volume Madry and Crumley (1990) present an 
example of GIS being used to explore the landscape settings of sites through viewshed 
analysis and the calculation of least cost pathways. The critique of GIS as encouraging 
retrogressive positivist approaches, a situation that was recognized and criticized early 
within the adoption process (Wheatley 1993; Gaffney / van Leusen 1995), led to a greater re- 
liance on visibility and movement studies as the main outcome of this humanist turn. Con- 
sequently both techniques have been developed considerably, with visibility studies refined 
to include cumulative, Higuchi and directional viewsheds (Wheatley / Gillings 2000), and 
movement studies similarly refined to include various sorts of surface approaches, not least 
accessibility (Llobera 2000). 

To use visibility or movement is now almost de rigueur in archaeological applications 
and it is the *or" that we wish to address here. Both approaches are seen as a way of hu- 
manizing digital models, avoiding the detached, “seeing everything from above” view of the 
world and situating the agent within the landscape. There is, of course, a strong whiff of 
technological determinism surrounding these techniques, as both are relatively easy to per- 
form as individual functions without any understanding or discussion of the underlying 
technical or theoretical assumptions. The theoretical shortfall here is that for human beings 
movement and visibility are not separate “procedures” or “functions” but are intimately 
linked to produce the world as we perceive it. 


* University of Oxford. 
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Here we draw heavily on the theoretical work of Gibson (1979) and Ingold (2000; 
2011). Gibson, a psychologist specializing in visual perception, argued that movement and 
visibility are interlinked and that we engage with landscape through a constantly changing 
series of vistas, his theory of “reversible occlusion.” He also introduced the concept of “af- 
fordance,” whereby landscapes depend on the interaction of perceiver and environment, 
which come together to offer opportunities for action. For example, gentle slopes afford 
movement whereas a sheer cliff does not. Ingold has added a phenomenological interpre- 
tation to Gibson's ecological approach to visual perception: he argues that visual perception 
never occurs from a fixed point but through a “continuous itinerary of movement,” so the 
“unfolding of a field of relations [is] established through immersion within an environ- 
mental and social context.” This also introduces the importance of the “social context” in 
addition to the environment: cultural markers within the landscape assume extra signifi- 
cance as influences on movement. In this paper we are inspired by, and try to model, sev- 
eral of Gibson's and Ingold's concepts based on the notion that movement is primarily vis- 
ual. This moves us away from the well rehearsed arguments of space versus place, where 
the concept of “place” is a culturally constituted but essentially static understanding of the 
world, and towards one based on pathways, wayfinding and life as a constantly changing 
network of cultural and ecological experiences that are lived through moving. 

Central to this approach, and to using GIS technology more widely, is an understand- 
ing of "scale," a surprisingly complex and slippery issue often taken for granted and not 
well theorized or often discussed in depth in GIS-based archaeology (Lock / Molyneaux 
2006). Here we recognize two different forms of scale (Dobres 2000), “analytical scale” 
and *phenomenological scale." Analytical scale is the realm of research, the recognition, 
measurement and analysis of patterns, including operational concepts such as ratio, frac- 
tion and spatial resolution. Phenomenological scale is the lived scale of being in the world 
and experiencing embodied interactions with the material world as part of everyday life. 
This difference reflects the quantitative/qualitative tension present in archaeological 
GIS and is something we try to address in this paper through the integration of method and 
theory. 


2 Previous work: hillforts of the ridgeway 


This paper builds on previous work that explored various methods of modelling visibility 
and movement along the Ridgeway, a prehistoric pathway in central England, and its rela- 
tionship with a series of Iron Age hillforts (Bell / Lock 2000; Lock / Pouncett 2010). The 
earlier paper combined two elements: the use of least cost pathways (LCPs) to determine 
the possible route of the Ridgeway in later prehistoric times, and viewshed analysis to de- 
termine whether the hillforts were located to be visible from the Ridgeway. The initial LCP 
calculated from a cost surface based solely on a function of slope — an unconstrained LCP — 
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“fell off” the Ridgeway following the initial path of least resistance. Because the Ridgeway 
runs east-west along a ridge top, it was decided to add a topographic bias which would more 
accurately model the human intentionality of walking the path, in effect attempting to 
model human behaviour rather than the behaviour of an algorithm. A function of aspect 
was added which prioritized east-west movement by decreasing the friction costs as the as- 
pect approaches 90° and 270? to the path. The resultant constrained LCP approximately fol- 
lowed the modern Ridgeway and formed the basis for further analysis. 

A series of binary viewsheds were generated at 250 m intervals along the Ridgeway 
with a maximum viewable distance of 3 km in order to model near and middle distance vis- 
ibility as someone walked along the Ridgeway. Viewsheds were accumulated to produce a 
visibility index and the location of hillforts was tested against the high visibility areas, in 
order to establish whether they were located to be visible, which, as it turned out, they 
aren't. This analysis highlighted the influence of topography for modelling movement and 
the importance of scale for modelling visibility. The second paper developed these themes 
further, introducing a multi-scalar approach which sought to better integrate modelling of 
movement and visibility through the concepts of corridors of intentionality, reflecting the 
long distance aim of moving from Point A to Point B, and intermediate waypoints, from 
Point A to Point B via Point C. These were underpinned by the concept of “affordance,” 
highlighting the relationship between the perceiver or active observer and the environ- 
ment, and were intended to shift the emphasis of modelling movement away from the “ef- 
ficiency" of accumulated costs towards the human intention that underlies movement. 

Slope is scale-dependent and, therefore, so too are any costs based on slope, and the 
LCPs which are generated on the basis of the accumulation of those costs. Landserf (Wood 
1996) was used to calculate slope for multiple scales of analysis," using window sizes to a 
maximum of 101 cells (equivalent to a planimetric distance of circa 1 km). Marked variation 
was observed in the LCPs for different scales of analysis, with changing correlations be- 
tween the LCPs and the Ridgeway in areas of varying topographic definition. Different to- 
pographic features were picked out at different scales of analysis, with some LCPs following 
the line ofthe Ridgeway closely and others picking out passes or cutting across the heads of 
valleys. Additional complexity was added to the movement algorithm to incorporate a func- 
tion of visibility based on cultural landmarks being visible as waypoints along the intended 
route. Corridors of intentionality were defined for each leg of the journey using directional 
binary viewsheds for a corresponding cultural landmark, in this case an intermediate hill- 
fort. “Cones” of vision were calculated by inverting the direction of movement and by re- 
classifying the directional viewshed based on the proximity to the hillfort. 


I  Landserf uses quadratic approximation to calculate surface parameters. A quadratic trend surface is fitted to 
the data points corresponding to the centroids of the cells that lie within an n x m window. Surface parameters are 
calculated for the data point corresponding to the cell at the center of the n x m window. 
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Two key themes emerge from the work on the hillforts of the Ridgeway: 1) scale de- 
pendency, i.e., that the perception of topography and landscape and the influence of top- 
ography and landscape are dependent on scale; and 2) uncertainty, i.e., that considerable 
convergence and divergence occur between individual “solutions,” so this is an exploratory 
process. Both of these themes are developed further in this paper using the concept of view- 
paths. 


3 Case study: the Danebury long barrows 


Our study area is within the 450 km? of chalk downland that has formed the focus for land- 
scape studies around the Iron Age hillfort of Danebury, close to the Wiltshire/Hampshire 
border in southern Britain. The pioneering work of Palmer (1984) provided the base maps 
incorporating his extensive transcription and analysis of aerial photographs, subsequently 
digitized and used in an early GIS-based analysis of the development of this prehistoric 
landscape including Neolithic, Bronze Age and Iron Age monuments (Lock / Harris 1996). 
Here we concentrate on a group of four long barrows in the north-western corner of the 
Danebury area (Fig. 1), verified by fieldwork as follows: 

— Barrow 13, Bulford 1 (Grinsell 1957, 139) 

— Barrow 14, Sheer Barrow (Crawford 1958, LB77) 

— Barrow r5, Milston 2 (Grinsell 1957, 142) 

— Barrow r7, Milston 39 (Grinsell 1957, 142) 

Long barrows are a distinctive class of Neolithic funerary monument, constructed dur- 
ing the mid- to late fourth millennium and early third millennium BC (Pouncett 2008). 
The visual characteristics of long barrows have provided a focus for GIS-based research, for 
example in order to explore territorial markers through the combination of binary 
viewsheds and Thiessen polygons (Lock / Harris 1996), and modelling visibility to and 
from individual barrows by binary and cumulative viewsheds (Exon et al. 2000; Wheatley 
1995). Here we concentrate on integrating visibility with movement between individual 
barrows at various scales. 

The dataset used here corresponds to an area 7 km by 7 km (centred on National Grid 
Reference SU 185 455) with the locations of the barrows extracted from Lock and Harris 
1996. The DTM (digital terrain model) is based on Ordnance Survey Landform Profile data 
with a ro m cell size and a vertical resolution of o.o1 m, obtained from Edina Digimap. 
A 5 km buffer has been applied to the dataset to minimize the impact of edge effects; this is 
more than double the mean distance to global horizon (see below) and assumes the maxi- 
mum range of effective visibility. The DTM was re-sampled to a 25 m cell size to allow rapid 
calculation of total viewsheds (see below). The topography of the area ranges between 68 m 
and 201 m OD (Ordnance Datum - above mean sea level) although the long barrows are all 
at the lower end of the range and occur in a variety of landscape contexts: 


GARY LOCK, MARIZA KORMANN, JOHN POUNCETT 


Key N 


; 285 mOD | 


Fig. 1 | The Danebury area showing the study area used here and the location of long barrows. The barrows used 


for analysis in this paper are shown in red and numbered. 


— Barrow 13 is situated on the lower slopes of a natural spur close to the apex of a tight 
bend in a river at a height of 96 m OD; 

— Barrow 14 is located just below the crest of a natural spur overlooking the eastern bank 
of a river at a height of 120 m OD; 

— Barrow 15 is located on the upper reaches of a natural spur along the western edge of a 
river valley at a height of 112 m OD; 

— Barrow 17 is situated on the northern edge of a pass between two river valleys at a 
height of 120 m OD. 

By modelling movement between these barrows we are not suggesting that their use 
was contemporary or that people did, indeed, move between them. The emphasis of this 
paper is on developing a methodology, although the locations of the long barrows are vari- 
able enough to raise interesting questions about the possible choices to be made if moving 
between them: is movement favoured along river valleys or along the adjacent upland 
ridges, or both under different circumstances? 
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4 An integrated approach: moving and seeing 


An important attempt to integrate visibility characteristics within the process of modelling 
movement is Lee and Stucky's work on “viewpaths” (1998), in which a viewshed con- 
tributes to the friction surface used to generate a least cost pathway. This concept simulates 
the influence of visibility on movement and is developed further below. The viewpaths gen- 
erated by Lee and Stucky were based on the calculation of two visibility surfaces or total 
viewsheds (viewgrids): 

1. Vg (General) — the number of cells which are visible from each cell (“looking out”); 

2. "Vy (Dominance) - the number of cells from which each cell is visible ("looking in"). 
These visibility surfaces or total viewsheds assume that the height of the observer and 

the height of the target are not equal, i.e., that the view "looking out" is different from the 

view "looking in" (see Fisher 1996, fig. 2). Four different types of viewpath were calculated 

using friction surfaces based on different characteristics of V and V: 

- Hidden Path - the path which is the least visible from other locations, with friction 
costs based on the lowest values of V; 

— Scenic Path - the path with the highest visibility, with friction costs based on the in- 
verted value of V¿; 

— Strategic Path — the path with the least possibility of being visible from other locations 
but maintaining the highest visibility from the path itself, with friction costs based on 
the minimum of the value of V, or the inverted value of V; 

- Withdrawn Path - the path that is both the least visible from other locations and 
has the lowest visibility, with friction costs based on the minimum ofthe values of Vp 
or Vo. 

Each of these four paths are tested against a benchmark Euclidean Path based solely on 
the horizontal distance and height difference between cells, so that comparisons of visibil- 
ity, invisibility and efficiency (length of path) can be computed. 


4.1 Multi-scalar analysis 


The Hidden and Scenic viewpaths proposed by Lee and Stucky (1998) are also imple- 
mented here using a modified form of total viewshed where visibility is restricted to a spec- 
ified scale of analysis or maximum visible distance (a technique elsewhere referred to as 
“fuzzy cumulative visibility analysis,”, “FCVA”, Kormann and Lock 2014). By restricting the 
scale of analysis it is possible to model the range of visibility that influences movement as 
an agent makes his/her way through a landscape, in digital terms moving from cell to cell 
across a DTM. The locations which can be seen change as the agent moves from one cell to 
the next and the dynamic nature of this change is reflected in the values of the visibility in- 
dices for the resultant viewshed. 
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Visibility indices are calculated for each cell in a DTM using a circular kernel with a 
radius equivalent to the maximum visible distance for the specified scale of analysis, so for 
example a 200 m scale of analysis is modelling a maximum distance of visibility of 200 m 
in any direction. The visibility indices are standardized by the maximum possible value of 
the visibility index for the specified scale of analysis in order to allow like-for-like compari- 
son between different scales of analysis (Fig. 2). Values of the standardized visibility indices 
range between o (none of the cells in the kernel are visible) and 1 (all of the cells in the ker- 
nel are visible), indicating areas of lower and higher visibility. Two viewsheds were calcu- 
lated for each scale of analysis, replicating the viewgrids which were used as friction sur- 
faces for the viewpaths proposed by Lee and Stucky (above). 

Rather than using arbitrary values for the scales of analysis, we here introduce the 
importance of skylines (or horizons) as a characteristic of visibility and, therefore, of 
movement. The importance of these, particularly with regard to astronomical alignments 
and territorial boundaries, has long been recognized in archaeology (Ruggles et al. 1993). 
Conventional GIS-generated viewsheds fail to take into account the visual impact of land- 
scape features relative to the position of skylines or horizons, i.e., points where the land 
touches the sky. Alternative forms of viewshed, including horizon viewsheds and local or 
global offset viewsheds, have been proposed to address this problem (Fisher 1991) al- 
though, because these are not standard push-button functions, they have not been widely 
adopted by archaeologists or other GIS practitioners. Here we use a variant ofthe horizon 
viewshed (Fisher 1996) to determine the minimum and maximum scales of analysis 
used below, with the minimum scale defined as the mean distance to the local (closest) 
horizon and the maximum scale defined as the mean distance to the global (furthest) hor- 
izon (Fig. 3). 

Euclidean Paths (cf. Lee / Stucky 1998 and above) were generated from Barrow 15 to 
each of the other long barrows in the study area (these are also used as the baseline below) 
using the Path Distance? and Cost Path3 tools in ArcGIS 10.0. These are based on slope dis- 
tance, taking into account horizontal distance and height difference only, and were calcu- 
lated using a constant friction surface (value = 1). Observer points were created at 100 m in- 
tervals along each of the Euclidean Paths for the purposes of calculating horizons. Lines of 
sight were generated for each of these observer points using z values equivalent to an ob- 
server height of 1.60 m and a target height of o.00 m. An angular increment of 1? was used, 
resulting in 360 lines of sight per observer point (36,720 lines of sight in total). The result- 
ant lines of sight were broken into visible and non-visible segments and the end points of 
the visible segments were defined as horizons. Distances to the horizons were aggregated 


2 http://help.arcgis.com/en/arcgisdesktop/10o.o/help/index.htmlt/ How. the. path. distance. tools. work/ 
009Z00000027000000/ (seen 29. 06. 2011) 
3 Mhttp://help.arcgis.com/en/arcgisdesktop/10o.o/help/index.html£/Creating. the. least cost path/ 
009Z00000021000000/ (seen 29. o6. 2011) 


VISIBILITY AND MOVEMENT: TOWARDS A GIS-BASED INTEGRATED APPROACH 


29 


Ez =e A T 


d 


SEMA 


"E 


Fig. 2 | Restricted viewsheds for the 7 km by 7 km study area (local = 200 m, global = 2 km and total = 
the whole area) compared against the binary viewshed for Barrow 15, red = o (low visibility / not visible), 
green = 1 (high visibility / visible), with north to the top. 
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Fig. 3 | The Euclidean Paths between the barrows showing observer points at roo m intervals from which local 
horizons (green) and global horizons (red) were generated. 
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for each observer point, calculating the minimum and maximum values, and these were 
averaged to calculate the mean distance to the local horizon^ (207 m, equivalent to 9 cells) 
and the mean distance to the global horizon’ (2,130 m, equivalent to 79 cells). Although the 
actual values determined here are specific to the study area, the method (Pouncett 2013) 
by which they were established has wider applicability. 


4.2 Multiple paths 


Having established appropriate scales of analysis, initial viewpaths were generated using 
the visibility surfaces (V¿ and Vp) proposed by Lee and Stucky (1998). Visibility surfaces 
were generated for viewpoints corresponding to the centroid of each cell in the DTM in- 
cluding the 5 km buffer (462,400 viewpoints in total) using the Viewshed tool in Arc- 
GIS 10.0.6 V, and V, viewgrids were obtained by inverting the observer (OFFSETA V, = 
1.6 mand V; = 0.0 m) and target heights (OFFSETB Vc = o.o mand V, = 1.6 m). The maxi- 
mum theoretical value for any cell in the total viewshed is 462,399, since a cell cannot be 
visible from itself, although in practice values range between 17 and 67,516 for V, and 65 
and 74,182 for Vp. This difference between the range of values for the two viewgrids high- 
lights the significance of reciprocity (Conolly / Lake 2006, 229-230). Furthermore, in both 
instances the range of values is significantly lower than the maximum theoretical value, 
which suggests that the 5 km buffer has successfully eliminated any edge effects. 

Perhaps unsurprisingly, the highest values in both the V, and V, viewgrids corre- 
spond to hilltops (peaks and ridges) and the lowest values correspond to valley bottoms 
(channels). The ridge running along the eastern edge of the study area, interrupted by the 
pass on which Barrow 17 is located, is the most highly visible topographic feature in the 
landscape. Whilst there are more areas of higher visibility for the V¿ viewgrid than for the 
Vp viewgrid, they are also more fragmented, with localized areas of lower visibility along 
prominent hilltops. More pronounced differentiation can be seen between the V< and Vp 
viewgrids at smaller scales of analysis; this is particularly noticeable at the local scale (more 
below). The V, and V, viewgrids were used to generate visibility-based friction surfaces for 
the purposes of calculating Hidden and Scenic Paths as defined above. 

To provide a comparison between the various methods outlined here and generate a 
series of images that will serve to stimulate discussion, we compare the Euclidean, Hidden 
and Scenic Paths with a more conventional LCP (based on slope) generated using the for- 


4 The end point of the first visible segment of a line of sight, also referred to as an intermediate horizon (Fisher 
1996). 

5 The end point of the last visible segment of a line of sight, also referred to as a skyline (Fisher 1996). 

6  http://help.arcgis.com/en/arcgisdesktop/10.o/help/index.html£/Performing visibility. analysis. with . 
Viewshed. and. Observer. Points/oo9zoooooov8000000/ (seen June 29. o6. 2011) 
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Fig. 4 | The Scenic and Hidden Paths between the barrows calculated using the total viewshed compared against 
the Slope and Euclidean Paths. 


mula? proposed by Bell and Lock (2000). Figure 4 illustrates the differences between these 
four paths from Barrow 15 to Barrows 13, 14 and 17, calculated using the total viewshed (as 
opposed to the restricted viewsheds based on the minimum and maximum scales of analy- 
sis, i.e., the mean distances to the local and global horizons respectively). 

Itcan be seen that there is a strong correlation between the Scenic and Euclidean Paths 
from Barrow 15 to Barrows 14 and 17, whereas the viewpaths from Barrow 15 to Barrow 13 


7  Cost- Tan(Slope?)/Tan(r)) 
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follow different routes, with the Scenic Path following the line of a ridge and the Euclidean 
Path cutting across the valley below. Likewise, there is a strong correlation between the 
Hidden and Slope Paths from Barrow 15 to Barrows 13 and r7, each following the line of a 
valley. In contrast the Hidden and Slope Paths from Barrow 15 to Barrow 14 follow different 
routes, one following the lines of ridges and cutting across valleys and the other following 
the lines of valleys and cutting across ridges. In both instances, the degree of correlation be- 
tween the Scenic and Euclidean Paths and the Hidden and Slope Paths would appear to be 
dependent upon whether the initial direction of movement is upslope or downslope. Pref- 
erence is given to higher ground in the case ofthe Scenic Paths to maintain a good view and 
to lower ground in the case of the Hidden Paths to stay out of view. 

Although the Euclidean, Slope, Hidden and Scenic Paths all represent the same intent, 
i.e., to move from one barrow to another, Figure 4 shows that there are multiple possibilities 
for the same outcome based on different factors. The Euclidean Path minimizes the surface 
distance and is, therefore, the most “direct” route regardless of slope, whereas the Slope 
Path minimizes the relative cost, usually thought of as “effort,” and so in effect is a longer 
but easier route. In contrast, the Scenic and Hidden Paths are influenced by visibility, the 
former maximizing “looking out,” which can be considered as “the view” and the latter mi- 
nimizing "looking in," so that by walking that route you are less likely to be seen. We are, of 
course, assuming here that people being seen is not important for Scenic Paths, but is so for 
Hidden Paths. Scale is inherent within these different paths and in the decisions made be- 
fore and during walking them. For the Euclidean Path the scale can be considered “global,” 
as it is based on knowing where the end point is and continuously minimizing the distance 
to it, whereas for the other three paths decision making is more “local,” as it incorporates in- 
formation based on slope or visibility at particular locations along the path. This can be re- 
garded as the difference between analytical scale and phenomenological scale highlighted 
in the introduction, the latter attempting to be a digital representation of perception. 

It is worth exploring scale a little further in terms of an agent located within a land- 
scape and how that agent perceives the near and distant landscape characteristics and the 
opportunities for movement. Previous work on the hillforts of the Ridgeway (Lock / Poun- 
cett 2010) has highlighted the importance of scale-dependency and how the value of slope 
is linked to the distance over which it is calculated, so that different values will be obtained 
by assessing the terrain closer to, or further away from, the observer. Since the friction costs 
used to calculate the Slope Path are based on those values, visibility and movement are in- 
trinsically linked. It has long been recognized that visibility is scale-dependent and here we 
expand on that insight through the use of two different scales of analysis. The paths shown 
in Figure 4 are based on the total viewshed, i.e., unrestricted vision across the whole area. 
In Figures 5 and 6 the Hidden and Scenic Paths are shown using the restricted viewsheds 
for scales of analysis based on the mean distances to the local and global horizons, i.e., 
maximum visible distances of circa 200 m and circa 2 km respectively. These paths are re- 
ferred to accordingly as Total, Local and Global Paths in the text below. 


GARY LOCK, MARIZA KORMANN, JOHN POUNCETT 


E 

Key 

€ Barrows 

—- Local 

— Global 
Total 

— Euclidean 


Fig. 5 | The Hidden Paths between the barrows calculated using the restricted viewsheds (local = 200 m and global 
= 2 km) and the total viewshed compared against the Euclidean Paths. 


Looking in detail at the Hidden Paths first (Fig. 5), we can see for the route from Barrow 15 
to Barrow 17 a close correlation between all paths except the Euclidean Path, with each path 
following the line of the same valley. The Global and Total Paths are almost identical, con- 
sistent with the "restricted" view from the valley bottom, whereas the Local Path hugs the 
base of the slope along the northern side of the valley rather than following the valley bot- 
tom. All of these routes offer less opportunity for being seen when compared to the alter- 
native route along the ridge-top taken by the Euclidean Path. Greater variation can be seen 
between the routes from Barrow 15 to Barrows 13 and 14, most evident in the case of the 
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Fig. 6 | The Scenic Paths between the barrows calculated using the restricted viewsheds (local = 200 m and global 
= 2 km) and the total viewshed compared against the Euclidean Paths. 
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Total Path from Barrow 15 to Barrow 13, which follows the line of the valley to the south of 
Barrow 15. This is in sharp contrast to the Global and Local Paths from Barrow 15 to Bar- 
row 13: the Global Path skirts the northern edge of the spur on which Barrow 15 is located, 
cutting directly across the top of the adjoining ridge, before following the line of a valley to 
the west of Barrow 15 to reach Barrow 13, while the Local Path follows a similar route but 
with subtle differences, skirting the northern edge of the spur before hugging the north- 
eastern edge of the adjoining ridge. The Global and Local Paths use more immediate topo- 
graphic features to provide cover from closer view, whereas the Total Path is concerned with 
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being seen from further away. Similar subtleties of difference can be seen between the 
routes from Barrow 15 to Barrow 14 at different scales of analysis. 

Perhaps not surprisingly the Scenic Paths, i.e., those maximizing looking out, produce 
equally rich opportunities for interpretation (Fig. 6). There is a good correlation between 
the Scenic Paths from Barrow 15 to Barrow 14, with all of the paths following a similar line 
along or across the same topographic features. Greater variation can be seen between the 
paths from Barrow 15 to Barrows 13 and 17. Again, the greatest discrepancies are with the 
paths from Barrow 15 to Barrow 13 where, although the Global Path corresponds closely to 
the Euclidean Path, the paths for restricted scales follow different topographic features. The 
Total Path follows the line of the ridge to the southwest of Barrow 15, while the Local Path 
follows the line of the valley below. Interestingly the Scenic Local Path and Total Hidden 
Path have very similar routes along the valley bottom, perhaps indicating that even for the 
latter the immediate topography is of most importance in restricting "looking in." It is also 
interesting that the Local Paths for the routes between all barrows deviate the most from the 
Euclidean Paths — an important point as it shows the influence of localized topographic fea- 
tures which are taken into account by the observer when moving. 


5 Discussion 


The analysis described above has shown that there is no single path between two lo- 
cations and multiple possibilities exist depending on the intentions of the agent: to gen- 
erate a single least cost pathway and claim significance for it as an “exact” route is a naive 
and reductionist attempt at modelling movement, especially when the influence of vis- 
ibility is taken into account. This is clear in Figures 4 and 5, which show a marked vari- 
ation in the degree of correlation between the Hidden and Scenic Paths from Barrow 15 to 
Barrows 13, 14 and 17. In this section we present some initial work that attempts to move 
beyond single paths and explores the idea of corridors of intentionality through probabil- 
ity surfaces, which show the differing affordance of movement and visibility across the 
landscape. 

In areas where there is a strong correlation between different paths (for example, 
Fig. 6, the Scenic Paths from Barrow 15 to Barrow 14), there is greater confidence that an 
individual cell will lie on the line of a particular path. Conversely, in areas where there is a 
poor correlation between different paths (for example, in Fig. 6, the Scenic Paths from Bar- 
row 15 to Barrow 13) there is lower confidence that an individual cell will lie on the line ofa 
particular path. Also important here is the notion of reciprocity, which is explicitly acknowl- 
edged by the distinction between the V, and V, viewgrids used as friction surfaces for the 
Hidden and Scenic Paths (see above). Depending upon the values ofthe observer and target 
heights used (here 1.60 m and 0.00 m), it cannot be assumed that Point A can be seen from 
Point B simply because Point B can be seen from Point A. Similarly, depending upon 
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whether an isotropic or anisotropic cost surface is used, it cannot be assumed that the least 
cost pathway between Points A and B will be the same as that between Points B and A. 

In the case of the Slope Paths, one possible solution would be to calculate a least cost 
corridor between each pair of long barrows rather than two separate least cost pathways. 
Least cost corridors can be calculated by adding the individual cost surfaces for two points 
together so that each cell in the grid has a value which represents the lowest combined ac- 
cumulated cost to get to both points. In the case of an isotropic cost surface, the least cost 
pathway will correspond to the cells within the least cost corridor with the lowest values and 
as such, the cell values can be inverted and standardized with values ranging between o and 
1 to reflect the probability that a cell lies on a least cost pathway. 

Here we extend the notion of a least cost corridor to generate probability surfaces that 
show the areas in which a least cost pathway or viewpath is most likely to exist. Least cost 
corridors were generated for the Scenic Paths from Barrow 15 to Barrows 13, 14 and 17 using 
the Corridor tool in ArcGIS 10.0.8 To construct the image shown in Figure 7, three least cost 
corridors were calculated for each pair of long barrows, one for each scale of analysis (local, 
global and total), and converted to probability surfaces by inverting and standardizing the 
raster datasets for each of the least cost corridors so that the values of the cells ranged be- 
tween 1 and o.? The individual probability surfaces were clipped to show cells with values 
=0.95 (i.e., cells with a 9596 or higher confidence level) and overlain using the Cell Statistics 
tool in ArcGIS 10.o* to create a single surface showing the maximum probability values for 
the Scenic Paths for the three combined scales of analysis (contrast Fig. 7 and 6). 

Marked variation can be seen in the resultant probability surface, with differential spa- 
tial patterning in the probability values along the line of individual paths. Several key points 
arise: all of the cells corresponding to the lines of Scenic Paths fall within the range of the 
clipped values for the probability surface, although probability values vary; the areas of hig- 
hest probability for the paths from Barrow 15 to Barrow 14, where the Scenic Paths are 
strongly correlated, are more extensive than those for the paths from Barrow 15 to Bar- 
rows 13 and 17, where the scenic paths are poorly correlated; the areas where the greatest 
difference between the paths for different scales of analysis or initial points of divergence 
occur correspond to cells which fall at the lower end of the range of probability values. 

It is clear from this example that the probability surface provides a much richer and 
more sophisticated understanding ofthe likely paths between the barrows than the individual 


8 http://help.arcgis.com/en/arcgisdesktop/10.o/help/index.html£/Creating a least cost corridor/ 
009Z00000024000000/ (seen 29. o6. 2011) 

9 Cells with a value close to 1 are more likely to lie within the least cost corridor for each pair of long barrows 
(i.e., cells which have a high probability of lying on a path between the long barrows). Conversely, cells with a value 
close to o are less likely to lie within the least cost corridor for each pair of long barrows (i.e., cells which have a low 
probability of lying on a path between the long barrows). 

Io http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/How_Cell_Statistics_works/ 
009Z00000088000000/ (seen 18. IO. 2011) 
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Fig. 7 | The probability surface for Scenic Paths at different scales, blue = low probability (lowest values equivalent 


to <=95%), red = high probability (100 %), with values clipped to show those of interest. 


least cost pathways or viewpaths. By looking at the correlation between the probability values 
and the underlying topography it is possible to explore the reasons for the convergence and di- 
vergence of paths between two locations at different scales of analysis. In terms of movement, 
visibility and perception of the landscape, this is an attempt to model the complex choice of di- 
rection on offer as an agent proceeds along a pathway to an intended destination. The choice 
will involve what can be seen in the immediate vicinity, in the middle distance and in the far 
distance, as well as the perceived topography and the affordance for movement offered at 
those different distance ranges, together with the known direction of the destination. 
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The methods of analysis which underpin the integrated approach to modelling move- 
ment and visibility advocated in this paper both compliment and extend existing ap- 
proaches to modelling movement or visibility. Probability surfaces, for example, are a more 
transparent representation of the results of modelling movement than least cost pathways, 
showing the areas through which a path is most likely to pass rather than a single solution 
which may or may not correspond to an actual path. Similarly, scale-sensitive methods em- 
phasize the importance of perception with regard to modelling visibility, which highlights 
the need to perform analysis at a scale which is appropriate to a particular landscape rather 
than at a scale which is entirely arbitrary. The analysis of the Danebury long barrows pres- 
ented above has demonstrated the interpretative potential of these methods, which inte- 
grate movement and visibility in order to model human behaviour in a more realistic 
manner and using the mean values of the distances to the local and global horizons to 
model that behaviour at a “phenomenological” rather than an “analytical” scale. 
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Path Network Modelling and Network of Aggregated settlements: 
a case study in Languedoc (Southeastern France) 


1 Introduction 


The possibilities and constraints on the movement of people and goods is a key factor that 
determines practices and structures space. The evolution of the communication network is 
therefore very important for understanding the evolution of a region. Our study aims to 
analyse the changes in the pattern of ancient aggregated settlements through communi- 
cation network modelling. The reconstruction of the ancient communication network is a 
means of understanding the evolution of this pattern not as a succession of creations and 
abandonments of settlements — point maps — but as a changing network. 

The study area is located in Languedoc (southeastern France) and extends from the 
Mediterranean coast to the hills ofthe hinterland (Fig. 1). In coastal regions, the communi- 
cation network is traditionally associated with sea and river freight. However, this paper 
deals with the road network at the regional scale, not with trade at Mediterranean level. In 
addition, it should be noted that the extent of navigability on local rivers in the periods 
under study is open to question; for example, the most important river of the study area — 
the Vidourle river — was navigable only during high water before its embankment from the 
plain from the Middle Ages onwards. The road networks have therefore played a major role 
in trade at the regional level in this area. 

Thanks to numerous previous studies, we have a good overview of the ancient settle- 
ment distribution and dynamic in most of the study area (Favory et al. 1985; Favory et al. 
1994a; Durand-Dastes et al. 1998; Nuninger 2002). However, knowledge of the ancient 
road network in this region is still deficient (Fig. 2). Research has long focused on the axis 
of long-distance communication connecting Italy to Spain through southern Gaul, known 
asthe Via Domitia (from the Rhone to the Pyrenees), also called the *Heraclean road" in the 
context of the prehistoric periods. In the study area, its itinerary is reported by ancient 
sources (literary and epigraphic) and most of its course has been reconstructed from its 
milestones (Castellvi et al. 1997). However, inter-regional roads — the “medium” level of the 
communication network — are still poorly known. Following G. Charvet, it is generally con- 
sidered that the major communication routes — linking Nimes to other capital cities — fol- 
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Fig. 1 | Study area and subset for settlement pattern analysis. 
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low approximately the course of present-day roads (Charvet 1873; Clément 2003; Monteil 
1999; Provost et al. 1999; Fiches 2002, 21, fig. 4). It is also assumed that transhumance 
paths and medieval salt and pilgrimage roads were used from Roman or prehistoric times 
(Clément 2003; Favory et al. 1985; Favory et al. 1994b; Fiches 2002, 80-81). At the same 
time, in some areas, map and photointerpretation studies carried out at local scale propose 
numerous possible Roman paths (Parodi et al. 1987; Bonnaud / Raynaud 1994; Raynaud 
2002), since the communication network seems to be well preserved in the landscape of 
this region: the archaeological sites are indeed often close to present paths. It is never- 
theless very difficult to define the status of these hypothetical segments of ancient roads in 
the ancient communication network: were they part of the regional network or local paths 
between farms? 

These reconstructions of ancient routes are based on regressive methodology; which 
starts from present-day elements of the landscape and selects tracks identified as being re- 


1 This kind of methodology uses the position of places of interest (river crossings, settlements), landscape 
analysis and toponymy. It results in the assembly of various segments of linear elements in the present landscape, 
interpreted as remnants of old routes. 
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lated to the periods under study. However, there is very little evidence to date the use of 
these old roads. Moreover, studies on roads have shown that the routes change, mainly be- 
cause of attraction by later settlements, and may undergo minor or major course changes 
(Vion 1989). Roads have a history and their course results from a long and complex evol- 
ution that combines abandonments, changes in status and reactivations (Vion 1989, 89). 
Thatis the reason why the number of proposed routes for a given regional road increases in 
areas where precise studies have been carried out (see box in Fig. 2). 

The knowledge yielded by a regressive approach is too scant to handle ancient road net- 
works: the chronology and status of known roads are uncertain and many major connec- 
tions remain to be located. Therefore, we adopted a constructive approach. Optimal path 
modelling, that is, simulating the connections between contemporary archaeological sites, 
helps to apprehend the whole communication network in its relative chronology and hier- 
archy. 

In contrast to the methods that have previously been used to study the road network 
in Languedoc, the aim of modelling is not to find the ancient roads in a precise manner, but 
to locate the communication channels that potentially result from the distribution of settle- 
ments. Optimal path modelling helps determine the most plausible path between two 
places, taking into account factors affecting mobility. In a further step, the computation of 
optimal paths between a set of places — optimal path network modelling — gives the struc- 
ture of a network; the multiplication of modelled paths highlights the communication 
channels which, according to the factors considered, are the most likely: that is, develop- 
ment in the form of roads is most likely where the intensity of trade is potentially highest. 

Our goal was not to study the numerous paths, built or not, that connect all settle- 
ments, but rather to trace the main communication channels that cross the region. This 
scale of analysis restricted the study to a specific level ofthe network and guided the choice 
of sites for path modelling. The regional level of the communication network has been 
modelled from the distribution of the agglomerations. Our selection of the major sites is 
based upon the collective research programme “L'habitat groupé gallo-romain et les ag- 
glomérations secondaires en Languedoc-Roussillon" (Fiches 2002). These sites can be very 
large or quite small (2 ha or less), nevertheless they all present evidence of specific status, 
such as the organizational structure of buildings, the presence of religious activity, epi- 
graphy, etc. It is acknowledged that agglomerations played an important role in the spatial 
organization of the territory of Nimes between the capital city and rural settlements (Fiches 
2002, 15). So occasionally we will call these sites “centres,” referring to centres of gravity of 
the economy or population. 
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Fig. 2 | Roads traditionally considered as Roman. 
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2 Optimal path implementation 


The factors and parameters taken into account in this work are mainly chosen according 
to the environmental characteristics of the interior of Languedoc: they are based on top- 
ography. The routes sought correspond to the least constraint path; in addition, the impact 
of a perceptual factor was tested: the field of view. To do so, two least cost path modelling 
methods were used. 


2.1 Optimal path factors 


From the perspective of travel by land, several categories of factors can be considered. First, 
we assume that the paths must be cost effective: they should enable fast movement, while 
energy consumption when travelling should be low. In fact, we assume that the difficulty of 
a journey was more important than travel time for the periods considered — especially since 
in the region the distance between the agglomerations in Roman times is only 25 km on 
average. 

The basic constraints on movement are considered usually as topography (slope), as 
well as the existence of barriers such as rivers, or surface conditions (vegetation, a built 
road, swamp, etc.). However, apart from the difficulties related to our poor knowledge of 
hydrography and surface conditions in ancient times,? we can consider these factors as 
secondary, to the extent that various construction techniques were used (bridges, land sta- 
bilization, forest clearance etc.). We do not aim to simulate the effect of movement of an in- 
dividual in a given context, but to assess the strong constraints imposed by a space travelled 
frequently by communities. From this perspective, we can consider that the topography is 
the primary element. 

Energy consumption also depends on the means of transport considered. To simplify 
the implementation, only pedestrian locomotion has been taken into account to set the 
least cost path modelling. Of course, the use of a cost function for wheeled vehicles would 
change the structure of the communication networks produced by theoretical modelling, 
due to different critical slopes (Herzog 2013). Nevertheless, we assume that a model based 
on pedestrian movement can already give a good assessment of the potential for the organi- 
zation of the communication network. 

Second, energy consumption is not the only variable that influences human travel 
(Llobera 2000; Wheatley / Gillings 2002, 151, 155). For this study, we have considered that a 


2 We know that in the study area river fords have moved significantly since the periods studied (Clément 2003, 
26), especially as a result of the hydrological system of this north Mediterranean karstic region. In addition, our 
knowledge of the palaeoenvironment is limited and it is still not possible to map ancient vegetation and wetlands 
reliably. 
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path providing a good field of view will be preferred by travellers transporting goods — inert, 
or live, such as herds — for orientation and safety reasons. We assume that it is one of the 
factors in the establishment of roads in elevated places (as better perceptual conditions 
would counterbalance the topographical difficulties). When implemented, the *visual" fac- 
tor is usually based on viewshed analysis (Lee / Stucky 1998; Lu et al. 2008). We however 
prefer to apprehend this factor through the field of view (FOV). Like energy consumption, 
this factor is conditioned by topography (see infra). 


2.2 Slope-dependent least cost path calculation 


Various solutions have been proposed to weight the distances depending on slope (for a 
comprehensive review of the various algorithms used in archaeology see: Van Leusen 
1999; Van Leusen 2002, 6.5-6.7; Wheatley / Gillings 2002, 154-156). Like others pre- 
viously (Llobera 2000; Podobnikar et al. 2004), we relied on biometric measurements: the 
energy consumption values relating to movement on slopes were based on the results of a 
physiological study made on mountaineers (Minetti et al. 2002). The range of slope gradi- 
ents given in Minetti's paper is quite broad, so we modelled intermediate cost values using 
a polynomial function of 3rd degree. 

The relative energy consumption calculated according to the slope value is an anisot- 
ropic parameter. Therefore, the slopes should not be considered as a set of fixed values, 
since the friction of movement is always dependent of the direction of movement: there are 
strong differences when walking up, down or along the slope. The problem of the differ- 
ence between energy consumption walking up or down the slope is sometimes treated par- 
tially by introducing an offset on symmetric functions, so that the cost of walking down the 
slope is smaller than walking on flat terrain (Verhagen et al. 1999; Van Leusen 2002, 
16.10-16.18). Some circumvent the issue by considering that a path is a round trip (De Silva 
/ Pizziolo 2001). However, it can be argued that, for a round trip, one would choose the best 
option by considering separately the costs for each direction; these costs cannot be reduced 
to an average cost, as has been proposed in some studies (De Silva / Pizziolo 2001). More- 
over, it is also possible to follow different paths depending on the direction of movement. 

Taking into account the direction of motion in path modelling has long posed techni- 
cal difficulties. It is not enough to offer different cost values for walking up and down 
slopes, it is still necessary to know the direction in which the slope will be traversed during 
each trip. Some commercial GIS software allows anisotropic calculations, such as the 
ESRI's Path Distance method, which can calculate the effective slope for each combination 
of two pixels from their respective elevation values ("vertical relative moving angle": ESRI 
2009). This effective slope is then transformed using a “vertical factor," which can be de- 
fined by a table correlating cost values with slope categories (in our case based from the 
energy expenditure values published by Minetti et al.). However, the integration of one or 
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Fig. 3 | Proposed procedure for ancient path modelling. 


more additional factor(s) affecting movement can be problematic, because this tool pro- 
vides only one way to combine the effects of the slope and any other cost component(s), 
namely multiplication. Theoretically, there is no dependency between the parameters we 
have chosen; the FOV (field of view) factor is not more important on a steep slope than on 
flat terrain. Therefore, addition is a more suitable method of combination than multipli- 
cation. 

Our approach proposes a specific anisotropic procedure: it considers the (approxi- 
mate) direction of movement during the calculation of the friction surface. This procedure 
consists of two main steps (Fig. 3), in order to predict the approximate direction of move- 
ment between two places which might differ significantly from a straight line. The first step 
consists in computing a temporary accumulated cost surface (cost surface 1). In the second 
step the approximate direction of movement is estimated by the calculation of the aspect of 
the temporary accumulated cost surface; this allows us to calculate the effective slope from 
which we performed a final accumulated cost surface using an isotropic algorithm (ESRI's 
Cost Distance), and then the least cost path. 

The calculation of the effective slope is based on the cosine of the difference between 
the aspect of the relief and the direction of the movement (aspect of cost surface 1). When the 
difference between the aspect and the direction of movement is equal to 90° (i.e. when 


PATH NETWORK MODELLING AND NETWORK OF AGGREGATED SETTLEMENTS 


49 


one isotropic method (costdistance) 


Altitude 


(m) 


N 
e 


4 8 12 16 20 24 28 km 
Distance (about 23,1 km) 


Esri's anisotropic method (pathdistance) 
00 


Altitude n 


(m) 


N 
e 


4 8 12 16 20 24 28 km 
Distance (about 22,7 km) 


Proposed least-cost path procedure 


200 
13 

2. 

<£ 

- 20 

CLE 4 8 1216 20 24 28km 
15 30 45 60 75 90 105 120 135 150 165 195 210 225 240 255 270 285 300 Distance (about 26,4 km) 

— Esri's isotropic method (costdistance) — Esri's anisotropic method (pathdistance) 

— + Idrisi 3.2 isotropic method (cost) — Idrisi 3.2 anisotropic method (varcost) 
— Proposed least-cost path procedure 


E. Fovet 2011 


Fig. 4 | Comparison of the different least-cost path procedures. 


walking along the slope) the recalculated slope is null, when this difference is equal to o* 
(i.e. when walking up the slope) the value remains unchanged, and when this difference is 
larger than 90° (i.e. when walking down the slope) the recalculated slope has a negative value. 


slope... = cos (aspect, direction) - slope (Zakšek et al. 2008) 


recalc 
The resulting paths are quite different from those obtained with the ESRI's Path Distance 
method (e.g. Fig. 4). They are “smoother,” exhibiting fewer topographical variations than 
the ESRI's modules. In the region studied, consequently, the proposed procedure tends to 
generate elevated paths, following the ridges. In contrast, the Path Distance method gen- 
erally produces quite straight paths, hardly influenced by the topography; this favours the 
transition into lowlands and narrow valleys. The analysis of all the modelled tracks (be- 
tween a set of agglomerations) showed that the communication channels produced by 
these two procedures represent different types of paths that we also find in the actual land- 
scape. 

We consider that the differences with the ESRI's anisotropic method do not under- 
mine our least cost paths procedure, although it is true that in very complex topographical 
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Fig. 5 | FOV as the average angle of horizon (drawing by K. Zaksek). 


conditions (with craggy relief) our procedure causes unnecessary detours. In these rare 
cases, however, the introduction of the field of view (FOV) factor improves the model sig- 
nificantly. 


2.3 Introduction of the field of view factor 


The field of view (FOV) has been defined as the average angle of horizon (Fig. 5). It was 
computed by searching the horizon at the maximum distance of 7 km in 32 directions on 
360°, with a methodology similar to skyview factor computation (Yokohama et al. 2002; 
Zak&ek et al. 2011; Kokalj et al. 2011). Areas with the average angle of horizon close to zero 
do not have many obstacles in their horizon (such as plains, ridges etc.), and areas with 
a high average angle of horizon provide a low field of view (as in narrow valleys or pied- 
monts). 

In order to obtain optimal paths that take into account both the effect of the slope and 
the impact of the field of view, the friction surface was calculated as a sum of the FOV and 
the relative energy consumption weights (see supra). However, the impact of the FOV on 
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human decisions about movement is not measurable. It is hence essential to evaluate the 
combination of the FOV and energy consumption factors. We computed several models, 
changing the weight of the field of view: 

In the first model, the FOV was set to 1; this model yields a very poor influence by the 
FOV, as its values are much lower than the energy consumption ones — FOV values range 
from 1 to 27.7, with a mean of 2.87, while those of energy consumption range from 24 to 
356, with a mean of 33. A second model was produced using a FOV weight of 5; the values of 
this layer are still two times lower than the energy consumption ones. In the third model, 
FOV weight was set to 10, which yields a strong influence by the field of view since the 
values of the two factors are then similar. In the fourth model FOV weight was set to 20; the 
values of this layer are two times higher than the energy consumption ones. 

The introduction of the FOV to the modelled tracks gathers the paths either on the 
ridges or in the centres of plains and basins. 


3 Assessment of the models 
3.1 Comparison with roads previously proposed as ancient 


Path modelling validation is a difficult issue. The relationship between network models and 
roads considered as ancient (e.g., presumed Roman) is ambiguous: if a modelled channel 
does not fit the presumed path for a given route, one cannot really know if that is because 
the model is erroneous or because the route traditionally proposed corresponds to another 
period or another level in the road network (local roads and inter-regional roads vs. regional 
communication channels). The status of roads implies specific routes, because the main 
purposes are different (e.g., supra-local roads do not aim to connect all local places), the 
level of engineering and the organization are variable and the modalities of transport can 
also be different (pack animals vs. wheeled vehicles for bulky goods). Inevitably, several al- 
ternative routes existed and our modelled networks are designed to locate only one of them. 
Thus, at the present state of studies on ancient roads in Languedoc, optimal path modelling 
cannot really be tested against the roads considered to be ancient; conversely, however, 
models can provide an analysis of the roads considered to be ancient. 

The Via Domitia is an interesting example for the problem of the confrontation be- 
tween the road network models and the known ancient roads. We know that this interprov- 
incial road was constructed during the 2nd century BC (see supra); however, it is now gen- 
erally considered that this was even then an adaptation of a pre-existing route (Py 1990, 
623; Clément / Peyre 1991, 15-25; Castellvi et al. 1997, 16; Py/Vignaud 1998). In our region, 
its path can be divided into two parts (see Fig. 2): the section Sextantio-Ambrussum and the 
section Ambrussum-Nîmes. Close to Nîmes, its precise location is subject to discussion 
(Fiches 1985, 136—138; Fiches 1997, 63; Laforgue et al. 1997, 24). 
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Fig. 6 | Length of the longest optimal path segments of the network that follow the Via Domitia 
(with a tolerance of 500 m). 


The Sextantio-Ambrussum section of the Via Domitia is followed, by and large, by the various 
optimal paths between these two towns (Fig. 6), with the exception of those from the model 
strongly influenced by our visual factor (FOV set at 20). In contrast, to model the path of the 
Ambrussum-Nimes section, it seems necessary to integrate the FOV with a significant 
weight — directing the optimal paths in the open plain. For the other models the only way to 
reconstruct the Via Domitia correctly is to take into account the small town of La Conda- 
mine as a stopping point. Archaeologists usually consider that the section connecting Am- 
brussum to Nimes has changed over time. According to J.-L. Fiches (Fiches 1985, 136-138; 
Fiches 1997, 63) its passage in the plain could not have been put in place before the layout 
of the cadastre “Nimes A” to which the road fits securely, that is to say from the rst century 
BC on. It was proposed as early as 1930 (J. Igolen) that the pre-Roman path passed via the 
oppidum of Nages, 5 km to the north ofthe Via Domitia (see Fig. 2, road A-1). This proposed 
pre-Roman route is more readily reproduced by the paths from our least cost procedure 
(without the integration of the FOV): optimal path segments are a bit longer, but mostly 
much denser. 

Concerning the rest of the road network, we observe that, in general, optimal paths fol- 
low approximately the tracks of the roads considered to be ancient} (Fig. 7), although they 
do not always correspond to dense channels. The communication channels produced by 
the Path Distance method correspond better to the presumed inter-regional Roman roads 


3 Only roads E, N, O and B are not well represented. 
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Fig. 7 | Length of the longest optimal path segments of the network which follow the roads 
previously proposed as ancient (with a tolerance of 500 m). 
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(roads C, D and F). Regarding the other types of road, some are better represented by this 
model (B-1, K, L and particularly the ancient salt roads M), while others are better repre- 
sented in the networks modelled using our least cost path procedure (G and particularly H 
and the transhumance paths I). Considering the specificities of the paths modelled by 
our procedure, it is not surprising that they are the only ones that fit to the ridge paths (H, I 
and La). 

These comparisons let us see that no model alone can approach the complex ancient 
road network; the use of several models completing each other is needed. In some cases 
historical reconstitution of a single road requires that different parameters be modelled. 
For example, the section Nîmes-Nages of the inter-regional road B is best fitted by our 
procedure when the visual factor is given no, or a very small, weight (FOV set at o or 1). Its 
Nages-Villevieille section is followed only by the Path Distance model. For the section 
Villevieille-Puech des Mourgues, models using our procedure are the best (especially with 
FOV set at 1 or 5). 

As a further step, the analysis of optimal path network modelling at regional scale 
shows that some inter-regional roads are likely formed by connecting up separate re- 
gional routes. Let us take the example of the road D, which connects Nimes to Anderitum 
(nowadays Javols in Lozére): its section between Nimes and the massif of Bois des Lens 
(in the north of La Jouffe) is represented in the majority of the models (Path Distance, 
models with FOV set at 5, 10 or 20) as a portion of the route between Nimes and Mus. We 
can therefore propose that this section ofthe inter-regional road D corresponds primarily 
to the regional roads from Mus to Nimes, to which the route from Anderitum would 
connect. 

The road C also seems to be composed ofthe regional routes from Nimes to Prouvessa 
and La Jouffe on one hand, and from Mus to Mauressip, Nage and Espeyran on the other 
hand. The connection between these two sections is poorly represented by modelling path- 
ways between the regional centres. It is therefore interesting to note that in the area where 
these regional networks connect (south of Prouvessa) this road makes a detour (Fig. 8). 
This area is also marked by a Late Roman funerary area and the proximity of a Roman and 
Late Roman settlement. This detour is obviously not caused by the topography and could 
well be an artefact of the evolution ofthe road network. In fact, the segment oriented south- 
east to northwest of this anomaly corresponds to the channels connecting the oppida La 
Jouffe and Mauressip, which serve a set of iron ore processing sites and high quality stone 
quarries in the massif of Bois des Lens (Fig. 8). 

We can therefore assume that the linkage of the two old roads has been influenced by 
the existence of other routes passing nearby. At present, there is no element that would 
allow us to fully understand this phenomenon of interaction between the different levels 
of roads (local, regional ...): only more data for the southwest to northeast communications 
would be able to support this hypothesis. In any case, the comparison between this road 
and the network models suggests that the axes of passages in this sector are complex; the 
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tracks could have swung between the attraction of long-distance networks and the logic of 
regional or local communications, and could well have had a story much more dynamic 
than is indicated by the regression analysis of the road network. 

The confrontation between the road network models and the roads considered as 
ancient raises interesting questions but, as we said above, it does not help to validate or in- 
validate the different models. 


3.2 Comparison with the settlement pattern 


Another way to evaluate the different regional road network models is to observe their con- 
sistency with the entire settlement pattern at different periods. Indeed, it can be assumed 
that the road network was a significant factor for the choice of location of many sites. For a 
farm, the vicinity of a major road made it possible to market its produce more easily and 
probably lowered the financial cost of transport. Latin agronomists paid special attention to 
the *commodity" of transport in the rural economy: 
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A handy road contributes much to the worth of land: first and most important, the actual 
presence of the owner, who will come and go more cheerfully if he does not have to dread 
discomfort on the journey; and secondly its convenience for bringing in and carrying out 
the necessaries — a factor which increases the value of stored crops and lessens the expense 
of bringing things in, as they are transported at lower cost to a place which may be reached 
without great effort; and it means a great deal, too, to get transportation at low cost if you 
make the trip with hired draught-animals, which is more expedient than looking after your 
own; furthermore, that the slaves who are to accompany the master will not be reluctant to 
begin the journey on foot. (Columella, RR I, transl. H.B. Ash 1941). 


In order to evaluate the models we calculated the frequency of archaeological sites* in a 
range of distance bands from the main communication channels, for four different periods 
which represent the major phases of evolution of the settlement pattern in this region: the 
5th century BC, the 2nd BC, the rst AD and the sth AD. To do this, we abstracted the main 
channels from the path networks obtained, using the density’ of paths. We started by map- 
ping the densest channels (which represent 5 or 6 paths, or 4 paths for certain periods), 
then we proceeded by completing gaps in these main pathways by following the most dense 
channels in between (usually 3 or 4 paths, or sometimes 2 paths). In certain cases, where 
two agglomerations were not connected with this procedure, we lowered the first threshold 
(3 or 4 paths). We restricted this analysis to the core of the study area (Fig. 1); this subset of 
452 km? is free of border effects and corresponds to the hilly zone ofthe interior of Langue- 
doc (the type of environment for which the model was developed) where the archaeological 
data are most reliable. 

Except for the sth century BC, a proportionally large quantity of sites is situated at less 
than 800 m from the main modelled channels, while in zones at more than 1200 m the 
number of archaeological sites declines (Fig. o). About 40% of the sites are situated in one 
of the first two ranges of distance in most periods. We note that this percentage is quite 
stable despite the strong variation in the numbers of sites between the different periods. In- 
deed, there are ten times more sites occupied in the rst century AD than during the 2nd 
century BC, but the proportions of sites inside the corridors are similar. 

For the sth century BC, the model which seems the most appropriate is the one built 
with our least cost path procedures without taking into account the visual factor (FOV 
weight set at o). For all other periods, the “Path Distance model” is the one that has most 
sites in the corridor of 400 m, while our least cost path procedures obtain better correlation 
between 400 m and 800 m when the visual factor is given no, or a very small, weight (FOV 
set at o or 1). Models that take into account the FOV significantly (weight of 5 or 10) are suit- 
able for the 2nd century BC and rst century AD, but less so for the sth century AD. 


4 Excluding those used as nodes to model the optimal path network. 
5 Defined by kernel density of optimal paths within a radius of 125 m. 
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Fig. 9 | Road network models and settlement pattern: proportions of sites located in distance 
bands from the main communication channels (plain colours) and proportions of landscape 
taken up by these zones (hatched colours). 
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A relationship appears between the types of sites? and the network models (Fig. ro). 
For the prehistoric periods, it is difficult to determine whether the various models concern 
specific types of sites, because the number of sites is low. One can nevertheless note that in 
the 2nd century BC ephemeral sites (class A) predominate in the models in which the vis- 
ual factor is given no, or a very small, weight (FOV set at o or 1), while class B is slightly over- 
represented in the vicinity of the * Path Distance model". 

In the 1st century AD, distributions of classes for each model are quite similar to the 
analysed corpus. One can nevertheless note that modest settlements (class A and B) tend to 
be over-represented in the models from our procedure, whereas this tendency is reversed 
for the “Path Distance corridors," where class E (higher status settlements) and C (medium 
size settlements occupied for a long time) are better represented. For the end ofthe 4th and 
during the sth century AD, distributions of classes for each model show stronger differ- 
ences. Higher-status settlements (class D and E) are better correlated with models without a 
visual factor, while smaller and short-lived sites (class A or B) are under-represented. Con- 
versely, in the vicinity of models that integrate the visual factor significantly (FOV multi- 
plied by 5 or 10), ephemeral sites (class A) are clearly over-represented while long-lived and 
large settlements (class C and D) are scarcer. 

Overall, these two sets of charts show that from the 2nd century BC on, we observed far 
more small and ephemeral sites in the vicinity of the network models that integrate the vis- 
ual factor and therefore favour elevated places or mid-plain, especially the model where the 
FOV is multiplied by 5. In contrast, the network of pathways that is relatively uninfluenced 
by the topography and which favours the transition into lowlands (produced with the Path 
Distance method) is more strongly related to large and sustainable settlements. This sug- 
gests that we modelled two types of network that seem to have had separate functions, since 
they not only reflect different ways of moving, but also appear to serve specific categories of 
archaeological sites. This phenomenon may be related to rural activities, such as pas- 
toralism, but the characteristics of settlements attached to these models should be studied 
in greater detail before proposing a firm link to this type of path network. 

These results let us consider the coexistence of two different kinds of network, each 
governed by a specific logic of movement (and probably with different aims and means of 
transport). Therefore, we decided to use the two most typical models to study the evolution 
of these regional path networks: a *model A" built with Path Distance and a *model B," 
which corresponds to the model that integrates the visual factor and which proved to be the 
most attractive, or specifically related, to ephemeral sites, i.e., the model that gives a mod- 
erate but significant weight to this factor (FOV weight set at 5). 


G This analysis uses the results ofa classification ofthe settlements (from the 7th century BC to the 7th century 
AD) made on a supra-regional scale (using a corpus localized in different micro-regions in southeast and central 
France) during the ArchaeDyn 2 programme (Bertoncello et al. 2012). This classification is mainly built on the 
highly discriminating variables “surface” and “duration of the occupation" of the settlements. 
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Fig. 10 | Road network models and settlement pattern: types of sites located in a corridor of 400 m from the main 
communication channels (hierarchical classification from ArchaeDyn2 programme). 
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4 Evolution ofthe optimal path networks 


The evolution from the 5th century BC to the 1st century AD of the distribution of the ag- 
glomerations produces only very slight changes in the configuration of the regional optimal 
path network in both models (Fig. 11 and Fig. 12). This apparent stability shown by the op- 
timal path models over these periods suggests that new agglomerations did not really 
change the structure ofthe road network but more likely completed it. The oppida of the sec- 
ond Iron Age (settled between the late 4th and early 3rd centuries BC) as well as the small 
towns created during early Roman times (from the end of the rst century BC and especially 
during the early 1st century AD) were established along the previous communication chan- 
nels (especially in model A). This suggests that the major settlements created during these 
periods were established in the vicinity of ancient roads. Thus these centres seem to have 
benefited from pre-established communication structures, and did not necessitate the con- 
struction of additional roads to trade and exchange with other important places. 

However, the establishment of these early Roman towns also seems to have influenced 
the appearance of (or at least reinforced the existing trend of) alternative routes to the pre- 
vious agglomerations. This tendency is more pronounced in model B. This phenomenon is 
particularly noticeable in the case ofthe oppidum of Ambrussum, around which the network 
model B forms four points of intersection within a radius of 3 km (Fig. 12). The existence of 
alternative pathways is indeed likely in this case, as it may be that travellers wanted to avoid 
the legal constraints (passage fee) probably associated with crossing the river by the bridge” 
at the foot of the oppidum. Therefore we can assume that alternative channels (type B) co- 
existed with tracks connecting the oppidum to the other regional centres (notably the Via 
Domitia linking Nîmes to Sextantio and Lattes). 

Beyond the overall stability of the network suggested by the optimal path models, the 
question that emerges concerns the role of new centres in the agglomeration network: were 
the new towns complementary to or competitors of the ancient ones? 

The cases of Ambrussum and Lunel-Viel are interesting in this respect. Lunel-Viel was 
established to the southwest of Ambrussum, in the middle of the 1st century AD, at the 
intersection of pre-existing communication channels of the model A (Fig. 11), and where 
the alternative channels in the model B meet (Fig. 12). This let us consider that the estab- 
lishment of Lunel-Viel resulted in channelling this second type of flow (B), probably on the 
track of road P (see Fig. 2 and Fig. 7). Therefore, Lunel-Viel would have had the effect of 
helping to better control the traffic circumventing Ambrussum. This hypothesis seems to be 
confirmed by the topography of this very small Roman town (Fig. 13): the road Cami 
roumieu (road P) penetrated from the east to the heart of the town and led to a large dead 
end furnished with monumental architecture, which suggests strong flow control on this 
road. 


7 On this Roman bridge see Vial 2003, 390-391, n° 1-1 
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Fig. 11 | Road network models from the 5th century BC to the 5th century AD: model A. 
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Fig. 12 | Road network models from the sth century BC to the sth century AD: model B. 
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Fig. 13 | Lunel-Viel: spatial extent of the Roman town and roads serving the site (extract from Ouriachi 2009, 
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From this point of view the new agglomeration seems to complete the role ofthe older 
centre. However, one may also wonder if, on the contrary, Lunel-Viel would not have rather 
had the effect of weakening the position of Ambrussum by reinforcing the alternative chan- 
nels, especially since Villevielle, which is located further upstream, developed strongly 
from the 1st century AD. 

On any scenario, there are parallels between the histories of Ambrussum and Lunel- 
Viel. Ambrussum witnesses great vitality in the 1st century BC, but the upper city was aban- 
doned around the beginning of 2nd century AD, and only the road station located at the foot 
of the hill, where the Via Domitia passed by, lasted longer. In contrast, Lunel-Viel was cre- 
ated in the middle of the rst century AD and witnesses a phase of intensification in the 
middle ofthe 2nd century AD (followed by stability until the middle ofthe 3rd century AD). 
Therefore, the decline of Ambrussum could be related to the potential competition exerted 
by Lunel-Viel (and other places?) at regional scale; or, at the least, the regional role of Am- 
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brussum in the road network could have been transferred to Lunel-Viel. The existence of a 
long-distance communication axis (the Via Domitia) does not seem to have been sufficient 
to maintain the economic vitality of the oppidum: the place turned into a stopping point 
only. Path modelling cannot explain the decline of Ambrussum, but shows the likely role of 
the road network in the relationship between this site and Lunel-Viel. This aspect was not 
fully taken into account before, because the question of the roads was only addressed by the 
analysis of the roads traditionally regarded as ancient, which recognizes no direct connec- 
tion between these two sites (see Fig. 2). 

This example shows that optimal path network modelling can give new perspectives 
on the changing pattern of agglomerations: it leads us to rethink the relationship between 
these major settlements, and to reassess their relative position in the trade networks. 

For the Late Roman period, the changes in the pattern of agglomerations do not seem 
to have a strong impact on the structure of the communication network either (Fig. 11 and 
Fig. 12). The abandonment or decline of several originally indigenous aggregated settle- 
ments between the end of the 1st century AD and the middle of the 21d century AD pro- 
duce only a decrease in the density of certain modelled channels, while the settling of St- 
Julian during the 3rd century AD and the reoccupation of Puech des Mourgues in the 
middle of the 4th century AD has only the effect of increasing the density of the east-west 
channels in the hinterland. This relative stability is partly due to the phenomena of slight 
displacements of population centres that do not affect the path network models — such as 
from Lunel-Viel to Lunel-Viel2, Lattes to Maguelonne, Le Castellas to La Cabanne and prob- 
ably from Espeyran to another nearby site not included in the models because of uncertain- 
ties (St-Gilles). In addition, if the site of Puech des Mourgues did not appear very well inte- 
grated in the previous regional communication networks, St. Julian is located near a very 
old communication channel, which is found in the various prior network models. This sta- 
bility suggests that abandonments or declines cannot be explained primarily by a restruc- 
turing of the regional road network. 

However, we must be cautious regarding the Late Roman period, which should be 
studied using a wider range of settlements to fully perceive the evolution of the network of 
exchange. Indeed, the end of monumental architecture and euergetism makes the major 
places ofthis period much less easy to apprehend from archaeological data than for the pre- 
vious periods, which are far better served by epigraphy and evidence of religious activity. 
The Late Antique central places are hence probably poorly represented in our modelling. 
Therefore, at this stage of the study, we will make only preliminary observations: 

We observe that the modelled channels for this period tend to pass via the location of 
former aggregated settlements that were not used in the calculation. The configuration of 
the general pattern of agglomeration would induce movement through (or close by) these 
sites anyway. The case of Ambrussum has already been discussed above: the road station is 
the only remaining settlement in the period studied, despite the fact that the place seems 
still to be well located in the regional path network; nevertheless its regional role may have 
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been transferred to Lunel-Viel. The case of Mauressip is very interesting insofar as it seems 
to illustrate a process of the survival of some roads near abandoned sites: when new settle- 
ments are established along a pre-existing network (Nages and Plaisance in this case), they 
continue to maintain the structure of the road network, even if the original places are no 
longer active. Note that in the model A the perpetuation of the east-west road near the 
former oppidum. (road L) seems due only to the existence of the Late Roman settlements 
(Puech des Mourgues and St-Julian); this raises the question of whether there was a process 
of reactivation of this track during the late Roman Empire. 


5 Conclusion 


The simulation of communication channels in this study was performed using several 
models. In addition to the use of common commercial GIS software, a procedure for cal- 
culating least cost paths has been proposed. In addition, the impact of a perceptual factor 
has been tested (the field of view). 

The analysis ofthese models — from comparison with the presumed Roman roads, and 
consistency with the entire settlement pattern at different times — led us to consider the co- 
existence of two types of networks: a network of pathways relatively uninfluenced by the 
topography and rather straight, which favours the transition into the lowlands, and a net- 
work which integrates the visual factor, which favours elevated places. This second type of 
modelled network seems to have a separate function from the first: it does not just repro- 
duce a different way of moving but also tends to serve a specific category of sites, namely 
those of very short duration. 

The two types of road network model can give new perspectives on the changing pat- 
tern of agglomerations. They show a slow construction of the road network from the Iron 
Age up to the early Roman period (1st century AD). Oppida of the second Iron Age and 
small towns created at the beginning of the Roman period fit into the networks modelled 
for previous periods: settled in the vicinity of previous communication channels, these 
agglomerations seem to benefit from pre-established communication structures. In this 
sense, the new major settlements did not really change the structure but completed it. 
Nevertheless, the models suggest that each introduction of new major settlements created 
or reinforced certain roads, or sections of roads, which offered alternatives to the old towns. 
Also, itis open to question whether these new agglomerations had the effect of channelling 
some ofthe alternative paths or weakening the position ofthe older towns in the communi- 
cation network. This approach puts a different perspective on old questions about the evo- 
lution of the network of regional centres, driven by opposing logics of integration and 
competition. 

For the late Roman Empire, however, the possible impact of the changes in the settle- 
ment pattern — characterized by the decline of several originally indigenous towns — on the 
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organization of the communication network are poorly identified by the modelling of rec- 
ognized towns. 

The confrontation between the road network models and the roads considered to be 
Roman has provided more questions than answers. At this stage of research on communi- 
cation networks, the modelled results give a glimpse into a phenomenon of integration be- 
tween different status of roads. The regional road network modelling shows that some 
inter-regional roads are likely formed by connecting up regional routes. In addition, some 
anomalies suggest that tracks could have swung between the attraction of long-distance 
networks and the logic of local communications, and could well have had a story much 
more dynamic than is indicated by the historical regression analysis of the road network. 
The intermediate level of the communication network (supra-local) — incorporating the 
role of large rural settlements — should be taken into account in order to test this hypoth- 
esis. 
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Finding Byzantine junctions with Steiner trees 


1 Introduction 


In archaeology, computational approaches to the study of movement practices have become 
more and more popular in recent years, especially when applied to prehistoric contexts 
(Bell / Lock 2000; Llobera 2000; De Silva / Pizziolo 2001; van Leusen 2002; Fábrega 
Álvarez / Parcero Oubiña 2007; Snead et al. 2009; Llobera et al. 2011). The main aim of 
these methods is not just to reconstruct paths and road systems, but to use these models to 
understand the establishment and transformation of a hierarchical network of settlement 
nodes, also including a focus on off-site land use practices. 

In this paper, we will introduce an example of the use of GIS-based modelling tech- 
niques for the study of the Roman-Byzantine communication system in the Roman Prov- 
ince of Cappadocia, in modern Turkey. We will argue that adopting a landscape archaeology 
approach to the analysis of an historic road system contributes to a deeper understanding 
of settlement dynamics (Bell et al. 2002). A GIS-based, spatial approach allows us to 
study settlements and pathways in a more formal manner and, by doing so, to address 
questions about ancient cognitive geography and past decision-making processes (Rock- 
man / Steele 2003). By systematically analysing the structure and phenomenology of the 
Roman provincial landscape, indicators of power and social order can be extracted from ag- 
ricultural centuriation and communication systems (Witcher 1998; Laurence 1999). In the 
cultural context of Roman and Byzantine communication systems, GIS methods are 
usually employed to reconstruct ancient routes in the landscape. These reconstructions are 
constrained by the possibilities of testing the model with material evidence of roads (re- 
mains, milestones, bridges; see GraRhoff / Mittenhuber 2009). A further limitation con- 
cerns the integration of path models created at a micro-regional scale with a macro-regional 
approach that considers the road system as a coherent planned network of connections. 

We therefore aim here to combine path models constructed on the basis of topo- 
graphic and environmental constraints at the micro-regional level — like slope, distance to 
natural resources, and visibility — with a model of macro-regional movement in a strategic 
network composed of nodes and junctions. In this network, movement is structured and af- 
fected by "cultural" constraints, like site density and distribution, and the position of stra- 
tegic administrative, economic, religious and/or military centres. 
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The integration of these two approaches is realized by combining a GIS approach for 
modelling segments of the system (paths), based on the environmental context, with a 
Steiner tree vector-based approach that takes into account the nodes and edges of the sys- 
tem. As an example, we use the documented inter-regional road system of Cappadocia that 
was established in the Roman period. It follows the previous pattern of macro- (Persian 
Royal Road)! and micro-regional (tribal villages) movement axes (Mitford 2000). Written 
sources and archaeological remains document this pattern of communication routes, con- 
sisting of strategic connections along a set of way stations. It can still be recognized in the 
pattern of Seljuk caravanserais in historical times (Hild 1977). 


2 The Roman and Byzantine road system in Cappadocia 


For this case study we focus on the methodological aspects of reconstructing and analysing 
the inter-regional road system of the Roman/Byzantine province of Cappadocia in Asia 
Minor (current Kapadokya in Turkey). The road system of Roman/Byzantine Cappadocia is 
an inter-regional communication system, from which we will analyse only one segment and 
we will focus on path modelling at the micro-regional scale. For this, we will integrate path 
modelling and network modelling approaches, and take into account topography-based 
constraints (Llobera 2000; 2001), as well as the structure of the macro-regional network. 

The study region (Fig. 1) is dominated by mountains and characterized by a continen- 
tal climate with snow from November to April in the area above the tree line. The area 
studied is a high plateau, enclosed to the south by the Taurus mountain range and cut trans- 
versely by the Antitaurus massif. The regional settlement pattern and the road system of 
Cappadocia are strongly structured and defined by topography. From the Bronze Age on- 
ward, the river plains were the main areas of settlement and acted as the primary axes of 
communication. Assyrians, Hittites, Phrygians, Lydians and Persians successively occu- 
pied the territory of Cappadocia, and were later followed by Seleucid, Roman, Byzantine, 
Seljuk and Ottoman rulers. The region has been subject to a centralized administration 
since the Persian period, and this clearly influenced the strategic construction and long- 
term use of the road systems. A strong relationship is evident, for example, between the 
Roman road system and the Seljuk caravanserais. 

During the Roman Empire, strategic routes were set up from the provincial capital of 
Caesarea (currently Kayseri; Greek: Kaisarea) to connect it to Melitene (currently Malatya), 
over a total distance of approximately 250 km.? Caesarea, the starting point, was a central 


1 Herodotus 5.52—54; French (1998). On the Cappadocian pathways cf. Herodotus 52.2.2, French (1998, 18) and 
Hild (1977, 99 and Abb. 68) on the route from Elbistan to Melitene. For a localization of the route on modern top- 
ography cf. French (1998, 28: Kayseri sheet: Kayseri-(Aydinlar)-Kamber-Akmescit). 

2 The modern road connection measures 339 km. 
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place from Assyrian times. It became a Hellenistic polis and then capital of the Roman 
province of Cappadocia and it was a central node of the road system in central Asia Minor 
during both Roman and Byzantine times. The connection we are studying was the south- 
ern alternative to reach Melitene, via Comana (Caesarea — Arasaxa — Comana — Kokousos — 
Melitene; Hild 1977, 85; Mitford 2000, 991). Along the route, way stations were built at dis- 
tances of approximately 40 km. A number of milestones, some of them found in situ, are at- 
tested along the route to Comana (Hild 1977, 85; French 1988 and 2012), and trace the cor- 
ridor used by the Roman road. The so-called Roman Itineraria list the distances between the 
stations along the roads. For this area, the Itinerarium Antonini (ItAnt) and the Tabula 
Peutingeriana (TabPeut) provide geographical information about the paths used. 

The pathway to Melitene passing by Comana seems to have been in use only until the 
Early Byzantine period (4th — 7th century AD; Hild 1977, 84). In the 7th century, an alter- 
native northern route passed via the fortress of Taranta to Melitene (for documentation and 
sources see Hild 1977, 84 n. 73). According to the sources, the area ofthe Antitaurus moun- 
tains, made desolate by the Byzantine campaigns against the Arabs, was re-organized dur- 
ing the 8th century. The eastern frontiers of the Byzantine Empire were restructured and 
the road system changed. In the Middle Byzantine period (9th — 11th century AD) a new al- 
ternative route diverged from Arasaxa - the first node of the Roman road and maintained as 
a Byzantine junction — to the north, passing via Tzamandos and Ariaratheia to Arabissos (as 
documented by Arab itinerary sources: Hild 1977, go and n. 108). This last connection is an 
alternative route that passes via Arabissos. Arasaxa still seems to have been used as a way 
station during the Mamluk campaign against Caesarea in 1277. The northern variant via 
Taranta was re-established in the Seljuk period and still serves today as the main route from 
Kayseri to Malatya (Hild 1977, 84, note 75). The southern variant is no longer used for mod- 
ern traffic (Hild 1977, 88). 


3 Case study 


For this paper, we only consider the first stretch of the route, from Caesarea to Arasaxa, the 
first way station on the route listed in the ItAnt. It was still in use as a caravanserai in Otto- 
man times. The straight-line distance between these sites is 34.87 km. In order to reach 
Arasaxa from Caesarea, the route has to move uphill and pass between the high Argaeos 
mountain to the south, and a lower mountain range to the north. On this stretch, the road 
has to cross several river valleys that are dry most of the year. As far as we know, no other 
settlements were connected by the road in Roman times.4 In Byzantine times, however, 


3 The extinct volcano of Erciyes Dağı, the highest mountain in central Anatolia, with its summit reaching 4013 m. 
4 The TabPeut links Caesarea first to Sinispora (?) (35.52 km) and then to Arasaxa (19.24 km; cf. Mitford 2000, 
990). The name Sinispora is corrupt, however (it mixes two other toponyms), and should be eliminated (Ramsey 
2010 [1890], 272). For an approach comparing Itinerary sources to reconstructed roads see French (1974). 
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Fig. 2 | Reconstruction of the stretch between Caesarea and Arasaxa (Hild 1977, map 6). 


a new military and religious settlement pattern developed along the road and probably rede- 
fined it. 

The reconstruction of the stretch by Hild (1977, 85 and map 6; Fig. 2) shows a route 
that passes through the sites of Moutalaske (currently Talas; a Byzantine 6th — 13th century 
hermitage site with religious buildings), Sari Han (a Seljuk caravanserai dating from the 
13th century, partly built from Byzantine spolia) and Meskiran Kalesi (currently Meskuan; 
a Byzantine castle). The sites of Iskokson (a small Byzantine site with a church) and Sakal- 
tutan (a possible Byzantine castle where a Roman milestone was found, not in situ) are 
indicated as additional waypoints (Hild 1977, map 6). About 4 km west of Sakaltutan, the 
route passes north ofthe excavated rock-cut monastic complex at Kepez (Hild / Restle 198r, 
map; Decker 2007, 242). Obviously, these are all post-Roman sites, and the reconstruction 
is therefore only pertinent for the Byzantine period.5 The sources attest to the changing pat- 
tern of communication routes in this region. The Byzantine and Medieval phases show a 
restructuring of the road from Caesarea to Arasaxa along a new settlement pattern consti- 
tuted by religious (Moutalaske is a birthplace of saints) and military landmarks (Sari Han is 
a caravanserai; Meskiran Kalesi a Medieval fortress). 


5 According to Hild / Restle 1981 the settlements in this area date from the 4th — 13th centuries AD. 
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The existing information about the diachronic settlement pattern, roads and com- 
munication system in this macro-region offers a convenient dataset for exploratory pur- 
poses. However, no comprehensive archaeological survey at the micro-regional scale has 
been undertaken across eastern Cappadocia.* Only some isolated surveys have been done, 
focusing on churches and Byzantine paintings, and almost all physical remains from the 
ancient world have been destroyed or re-used. Apart from the displaced Roman milestone 
at Sakaltutan, the only known traces of the Byzantine road have been discovered near Mes- 
kiran Kalesi (Hild 1977, 85). Modern travel journals from the r9th century (see Mitford 
2000, 985-986) provide some additional descriptions locating places and roads (Ramsey 
2010 [1890]). However, these are not always reliable. 

In order to analyse the course of the paths, and their stability or possible change over 
time in relation to the documented pattern of settlements, we have used the following da- 
taset including topography, historical geography, archaeological and epigraphic sources: 
- Digital elevation model: Aster DEM at 3o m resolution; 

-  GUGK (Glawnoe Uprawlenie Geodesii i Kartografi) Maps (scale 1 : 200,000) 

— Sites, bridges, road remains: 
° Google Earth-based localizations of ancient sites on the basis ofthe modern names; 
e Hild / Restle (1981); 

- Milestones: French (1988, 2012); 

— Road system description: Hild (1977); 

— General information: Map 64 Caesarea-Melitene (Mitford 2000) in the Barrington 

Atlas of the Greek and Roman World (Talbert 2000). 

— Itinerary Sources: 

e Itinerarium Antonini (Cuntz 1990); 

° Tabula Peutingeriana (Miller 1988). 

The ancient geographical dataset of the Roman Itineraria contains explicit information 
about settlements and distances as well as implicit information about the regional road net- 
work as a coherent infrastructure, connecting central places and a pattern of secondary 
settlements or stations. Archaeological information about the nature and typology of the 
secondary settlements related to the road system is not available for the region. Stations in- 
dicated in the Itinerary sources cannot always be located on the ground. 


4 Roman road building and least cost path modelling 


In the current case study we are dealing with the construction and maintenance of the 
roads that were administered by the Roman and Byzantine Empire. The archetypical 


6  Foran archaeological survey focusing on roads and paths in Cappadocia in the Sivas region, see Tuba Ókse 
(2007). 
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image of a Roman road is that of a straight route, paved with slabs of stone. The Roman 
administration indeed placed special emphasis on the efficiency of connections. The 
emperor Vespasian, for example, when advancing into Galilee in AD 67, records that he 
straightened and levelled existing tracks (Hucker 2009). All over the Roman Empire 
we find evidence that the construction of the major roads was aimed at creating speedy 
and reliable connections that could also be easily maintained. Even though constructing 
a Roman road must have been an expensive venture, once the roads were built they 
proved to be highly persistent, as is witnessed all around the former Roman Empire, 
where many Roman routes are still in use. However, in many instances the landscape 
offered challenges to Roman engineers, like steep slopes or wet areas that made the con- 
struction of straight stretches difficult (Quilici 1995). The Romans therefore adopted 
a flexible approach to road construction. Evidence for complex engineering solutions 
can be found in many places, in particular bridges, dikes and switchbacks, and some- 
times even tunnels. All in all, the Roman road system is more variable than is often 
thought. 

Where to put the roads must have been dictated by military, political and economic fac- 
tors, for which the speed, reliability and comprehensiveness ofthe connections between ad- 
ministrative centres must have been of primary importance. Two aspects of Roman road 
building make it somewhat different in terms of route (least cost path) modelling than most 
other cases considered in archaeology. These are connected to the primarily military func- 
tion ofthe roads. In fact, most major roads in the Roman provinces were initially built to ac- 
commodate the army advancing into enemy territory, and were only later fixed in place as 
“imperial highways." The roads had to accommodate not just the soldiers and horses, but 
also the supply trains of wagons and carts. 

From a modelling perspective, this poses a challenge: while extensive (experimental) 
research has been done to define equations that adequately model speed and energy expen- 
diture for travel on foot (see e.g. Herzog 2013), similar data are not available for wheeled 
transport using transport animals. One ofthe few authors who has made an effort to collect 
this information is Raepsaet (2002). He produced an equation that gives the traction force 
needed to get a cart moving, depending on the weight of the cart, the paving used, and 
the slope.? If we use this equation with the figures for a typical Roman cart, the carpentum, 


7 On the organization and infrastructure of the cursus publicus see Kolb (2000). 

8 For a short introduction to Roman land transport based on literary, epigraphic and legal sources see Meijer / 
van Nijf (1992, 136—140). 

9 T=kP+Pi 

where 

T = traction force needed for movement 

P = weight of loaded cart in kg 

k = rolling coefficient 

i- slope in m/m 

The rolling coefficient k is composed of a pavement friction factor and an axle friction factor. 
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Fig. 3 | Ancient road remains in Cappadocia (Hild 1977, Abb. 62: road remains west of Elbistan; 
Abb. 68: Byzantine bridge west of Akcadag). 
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drawn by 2 mules with a supposed maximum load of 500 kg (Roth 1998, 208-212), it in- 
dicates that it would be impossible to get these carts moving up a slope of more than 9%. 
While we can assume that in the case of steeper slopes wagonloads would be reduced, most 
Roman roads never take slopes > 15% (Hucker 2009; Quilici 1995). “Normal” least cost 
path calculations using hiking functions, like those defined by Tobler (1993) or Minetti et al. 
(2002) are therefore not best suited for Roman road modelling. 

Secondly, the positioning of Roman roads in various parts of the empire seems to in- 
dicate that visibility was an important consideration as well. In Britain for example, Roman 
roads consistently follow ridges and plateaus instead of valleys when given the choice 
(Hucker 2009). This may only partly have been for reasons of preventing ambush: when 
setting out the road, Roman engineers (gromatici) would have had to set up survey stations 
and sighting points on the route that should be intervisible. Typical distances between sur- 
vey points would be in the order of 2-3 km, but if braziers were used longer distances could 
be measured as well (Hucker 2009). The Roman road system in our study area was prob- 
ably constructed in the same way. 

Much less is known about the construction and maintenance of the Byzantine road 
system. The major waypoints were still maintained as traffic nodes, but instead of emphas- 
izing speedy connections between the major administrative centres, the Byzantine road 
system constitutes a more fine-grained transport network that also connected the minor re- 
ligious and military centres in the area. In eastern Cappadocia we have several archaeologi- 
cal remains that relate to the Byzantine road system (Hild 1977; Hild / Restle 1981): in the 
area of interest a bridge on the road from Kayseri to Malatya is documented (Hild 1977, 
Abb. 68). From the material remains it is also clear that difficult topography was sometimes 
preferred over wetlands (Hild 1977, 95 and Abb. 62, 68; Fig. 3). 
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5 Modelling the Roman route with least cost paths 


Because of the lack of applicable cost functions that take into account wheeled transport, we 
have compared the two most commonly used hiking functions (Tobler's and Minetti's) and 
a least cost function avoiding all slopes over 9% (which might have been prohibitive for 
wheeled transport). We have calculated least cost paths between Caesarea and Arasaxa 
based on these functions using a slope map derived from the Aster DEM.'? The modelled 
paths between Caesarea and Arasaxa all show similar routes, passing north of all the Byzan- 
tine sites mentioned as waypoints by Hild (Fig. 4). The distance of the Tobler and Minetti 
paths is 37.26 and 37.5 km respectively, and both routes are relatively straight, with a sinu- 
osity" of 1.07 and 1.08 respectively. According to the ItAnt (Cuntz 1990) and TabPeut (Mil- 
ler 1988), the length of the route from Caesarea to Arasaxa measures 35.52 km (24 Roman 
miles (mp); 1 mp = 1.48 km (French 1998, 146, fig. 2). The walking time needed to follow 
the shortest (Tobler) route is approximately 9 hours and 15 minutes,'? and would therefore 
fit well within a day's travel. The route avoiding all steep slopes is 39.69 km long (sinuosity 
1.14, approximately 10 hours walking). 

These figures are in line with the available information on travel speeds in the Roman 
Empire. For pedestrian movement, Kolb (2000, 310, tab. 1 and note 1; ibid. 311) gives an 
average of 20-25 mp (30-37 km) per day. Similar figures are supplied for transport using 
pack animals like mules, donkeys, horses and in some regions even camels and drom- 
edaries (on average 24 mp or 35 km per day; Kolb 2000, 312, tab. 2). Higher speeds could of 
course be achieved by riders on horseback or carts, especially with a regular change of ani- 
mals. 

The sources are less clear about travel speeds of heavy transport with wagons and carts. 
For ox-carts, a speed of 8 mp (12 km) per day is given by Kolb (2000, 316, Tab. 5). This, how- 
ever, seems to be on the low end of the scale. For the Early Medieval period Bachrach (1993, 
717) estimates a maximum speed of 9 mp (15 km) per day. Roth (1998, 211), however, 
reports that 1gth-century American ox-carts could cover 19 to 24 km (12 to 15 mp) per day. 
Mule-carts are faster and could easily cover 19 mp (30 km) per day (Bachrach 1993, 717). 
Since there are no equations that specify the relationships between slope and travel speeds 
for various kinds of animals and carts, it is impossible to be more specific about the effi- 
ciency of the modelled paths with regard to the various modes of transport. In either case, 
it seems impossible for ox-carts to cover the whole distance between Caesarea and Arasaxa 
in one day, and even for mule carts it seems improbable. 


10 The paths were calculated in ArcGIS 9.2, using the Path Distance module in order to account for the effect of 
anisotropic slope (see Herzog 2013). 

11 Thesinuosity index (Mueller 1968) is a measure of the directness of a route, it calculates the deviation from a 
straight line by dividing the route's length by the straight-line distance. 

12 The walking time calculations are based on Tobler's hiking equation. 
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Fig. 5 | The Derevenk valley (http: //static.panoramio.com/photos/original/20000398.jpg; 
copyright: Efkan Sinan). 


The least cost path models do not adequately account for the crossing of the Derevenk val- 
ley, which has very steep sides (Fig. 5). The Aster DEM allows a relatively easy crossing just 
east of the site of Moutalaske, which seems a highly unlikely solution given the steep 
descent. The vertical accuracy and horizontal resolution of the DEM is probably not good 
enough to reflect the actual topography of the Derevenk valley. Obviously, by taking the 
route further up- or downstream, the building of an expensive bridge may have been 
avoided. Forcing the least cost path model to avoid crossing the Derevenk valley creates a 
more southerly variant (Fig. 6). However, it still does not pass through any of the sites men- 
tioned by Hild. We can therefore conclude that Hild's reconstructed route, even when we 
force the least cost path to avoid the Derevenk valley, is not the most efficient option to reach 
Arasaxa from Caesarea. Hild's route has a length of 40.97 km (approx. 11.5 hours walking, 
sinuosity 1.17). 

This is also true when we take into account the Byzantine waypoints. Making the least 
cost path model pass through the Byzantine sites involves creating paths from site to site, 
rather than modelling a single route from Caesarea to Arasaxa. This results in a path with a 
total length of 39.17 km (approx. 1o hours walking, sinuosity 1.12). Avoiding the Derevenk 
valley takes 40.39 km (approx. 10.25 hours walking, sinuosity 1.16). Following Hild's route 
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from Caesarea to Arasaxa would therefore add at least one hour extra walking. But since his 
reconstruction is plotted on 1:800,000 scale topographic maps, it is probably not very ac- 
curate. 

By comparing different sources we can conclude that the use of least cost path 
methods for road reconstruction conflicts with the unpredictability of past decision-making 
processes as a result of non-optimal environmental and/or cultural choices. However, we 
can gain a better understanding of these processes using a model based on formal criteria 
that involve quantitative and qualitative aspects related to environment and settlement his- 
tory. 

For example, if we include visibility into the model, by calculating a total viewshed of 
the area and letting this progressively weigh as a cost factor, the modelling results indicate 
that the routes with best visibility are found further to the north (Fig. 7). This suggests that 
Roman road building in the area was probably not dictated by optimal visibility. It also 
prompts the question if there was something like “good enough" visibility for the Roman 
gromatici, and whether using a total viewshed would be the best way to model this (see also 
Verhagen / Jeneson 2012). 


6 Junctions and Steiner trees 


Hild's reconstruction shows a junction that is located approximately halfway between the 
sites of Phlabianai (a Byzantine church site, 4th — 13th century AD) and Iskokson. Here, a 
second route branches off to the southeast, leading to the Byzantine settlements of Cebir 
and Tomarza and onwards to the Late Roman/Early Byzantine site of Kiskisos (Hild / Restle 
1981, 206). 

Additionally to slope, which influences the feasibility and stability of the paths, a struc- 
turing role in movement is played by nodes and junctions (Gibson 2007). Junctions are im- 
portant elements in a coherent road system, but have been neglected in many case studies, 
in part because few applications of GIS in the field of Roman and Medieval archaeology are 
available (Witcher 1998; Bell et al. 2002; Fiz/Orengo 2008; Bellavia 2006; Gaffney 2006; 
de Soto / Carreras 2008; de Soto 2010). A GIS-based path modelling procedure using only 
archaeological sites as waypoints does not allow for the creation of junctions that optimize 
movement along the whole strategic network. From a mathematical point of view, the prob- 
lem of finding the optimal interconnection of a set of points is covered by the Steiner Tree 
Problem, a generalization of the Minimal Spanning Tree Problem. Given a set of nodes 
(points) with edges (connecting lines) that connect various pairs of points and have associ- 
ated costs, the minimal spanning tree seeks the least cost tree (i.e. a collection of edges hav- 
ing no loops, and the points they connect) containing each of the given points. In a Steiner 
Tree Problem, additional points (Steiner points) may be utilized if they can help lower the 
cost of the tree (Fig. 8). To avoid confusion, the original points that must be connected are 
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Fig. 7 | Least cost paths between Caesarea and Arasaxa, based on Tobler's hiking equation and different weights for visibility. 


T1 
S1 S2 
Fig. 8 | Steiner tree connecting points Tr, T2, T3, and T4 by using 
Steiner points Sr and S2. The four terminal points could be connected 
without using S1 and S2, but the cost (i.e., length ofthe tree in this case) 
would be higher. T3 


referred to as "terminal nodes." These problems are often situated in the Euclidean plane, 
in which case costs are the Euclidean distances. Since distances between settlements in real 
landscapes should not be measured as Euclidean distances, but as (non-uniform) cost dis- 
tances, an additional complication is found in creating the optimal network connections. 
The Euclidean Non-Uniform Steiner Tree Problem, described by Frommer / Golden 
(20072), aims to efficiently solve the Steiner Tree Problem for a network of points using 
non-uniform distances. In the Euclidean Non-Uniform Steiner Tree Problem, each location 
in the Euclidean plane has an associated cost. An edge connecting two points will have a 
cost that depends on the cost of the locations through which it passes (a least cost path). De- 
vising a network of paths connecting sites in a rugged mountainous landscape can be for- 
mulated as a Euclidean Non-Uniform Steiner Tree Problem, with location costs dependent 
on land use, slope, elevation and other factors. 

Along these lines, the problem of trying to find possible routes between archaeological 
sites given known waypoints can be formulated as a Euclidean Non-Uniform Steiner Tree 
Problem. The existing sites can be thought of as the terminal points in the tree that must be 
connected in a reconstruction of an ancient network of paths or roads. Steiner points help 
tolower the cost ofthe overall tree, i.e., they result in more efficient pathways. Furthermore, 
Steiner points always have at least three edges leading to/from them, and hence may func- 
tion as junctions in a (road) network. Generating Steiner trees for larger sets of points is 
computationally intensive, however. Furthermore, no implementations of algorithms that 
generate Steiner trees from point data are available in GIS.» 


13 The v.net.steiner module in GRASS (grass.osgeo.org) only finds the optimal connection between a subset of 
points in an existing network (i.e. with predefined edges); it will not create a new network. 
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7 A genetic algorithm to solve the Euclidean non-uniform Steiner tree problem 


Frommer and Golden (2007a) developed a genetic algorithm to solve the Euclidean Non- 
Uniform Steiner Tree Problem, and applied this to sample problems with structured land- 
scapes (e.g. hills, rings, etc.). The algorithm was later combined with GIS data in an exer- 
cise aimed at connecting recreational sites in Stowe, Vermont (Frommer / Golden 2007b). 
A weighted combination of layers of GIS data (elevation, land use and slope) was used to 
specify the underlying cost structure. 

The problem is situated on a hexagonal grid. Each hexagonal cell on the grid may con- 
tain at most one point, and has a cost associated with it. Given two points, the edge connect- 
ing them is defined as the shortest path between them, which is not necessarily the straight 
line segment connecting the points. The cost of the edge equals the sum of the costs of the 
cells through which the edge passes plus one half the costs of the two end points' cells. 
Using a grid reduces the size of the solution search space from infinite to finite, though it 
is still potentially very large. It also allows the problem to be represented in network 
(or weighted graph) form. A genetic algorithm is only one possible solution, and some of 
the trade-offs between the genetic algorithm and other approaches are discussed in 
Frommer / Golden 2007a. What is important for this work is that the genetic algorithm 
finds good (and in some cases optimal) solutions for reasonably sized problems relatively 
quickly. It may not necessarily be the fastest or best performing algorithm, but it is quite 
flexible to changes in problem formulation. 

Genetic algorithms are based on biological evolution and natural selection. They typi- 
cally consist of a population of individuals representing solutions that change over time 
through the application of crossover and mutation operators. A fitness function assigns a 
value to each individual in the solution with regard to how well it meets the problem goal. 
In the genetic algorithm we employ, each individual in the genetic algorithm population is 
defined by a set of potential Steiner points. The fitness function is the cost of the minimal 
cost tree connecting the terminal points and possibly using some of the Steiner points. 
Lower cost solutions are chosen with higher probability than higher cost solutions, and are 
then used to generate new individuals through the use of crossover and mutation. This gen- 
etic algorithm is discussed more fully in Frommer / Golden 2007a. 


8 Application 


The genetic algorithm was applied to find the minimal spanning tree that connects the way- 
points on the route from Caesarea to Arasaxa and those on the route branching off to the 
southeast towards Cebir via Phlabianai, as specified by Hild (1977). The results are shown 
in Fig. 9. The genetic algorithm solution does an effective job of connecting the terminal 
points, while avoiding high cost regions wherever possible. The deviations from the Tobler 
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least cost paths calculated between the way points are clear. The two detours to the north 
(right from the start in Caesarea, and between Iskokson and Sakaltutan) are more in line 
with Hild's reconstruction than with the least cost path. The junction near Phlabianai is 
close to where Hild hypothesized it. Note that the genetic algorithm finds a detour in the 
upper left thatis not found in the other models. This solution uses a narrow low-cost path of 
approach to Moutalaske that yields a cost approximately 3% less than the next best route. 
A question to consider is how likely people would be to find the less intuitive detour when it 
only provides relatively small savings. Such considerations can be incorporated in the algo- 
rithm, as described below. 

Since the algorithm works on a hexagonal grid, the number of directions of edges is 
limited to six. So solutions sometimes exhibit staggered (zigzag) paths when they need to 
move in a direction other than one of the six. Because of this, the efficiency of the calculated 
route appears to be very low. The connection created between Caesarea and Arasaxa has a 
length of 50.84 km (approx. 13.5 hours walking, sinuosity 1.46). The zigzags could be re- 
placed by the lines they are attempting to follow in a post-processing routine. Given that 
this work represents only an initial exploration, the post-processing has not been imple- 
mented. Alternatively, the current set-up could be replaced with a rectilinear grid and edges 
allowed to take on any direction. Cost would be calculated using numerical line integrals. 
Steiner point locations would still need to be restricted, possibly to the centre of each cell. In 
addition, the grid could have multiple resolutions, with a finer grid used in areas of highly 
variable cost structure. 


9 Concluding remarks 


We have shown that the traditional application of least cost path algorithms will only calcu- 
late optimal pathways between pairs of points, and are not effective for finding intermediate 
waypoints, like those hypothesized in Hild's reconstruction of the Byzantine route between 
Caesarea and Arasaxa. The only way in which least cost paths will approach the location of 
the supposed route is by progressively including known waypoints and calculating least 
cost paths between them. Even then, the resulting path network will not be the optimal sol- 
ution to connect a set of known points. It also deviates substantially from the route sug- 
gested by Hild. We applied the Euclidean Non-Uniform Steiner Tree Problem using the 
genetic algorithm developed by Frommer and Golden (20072, 2007b) to the issue of find- 
ing the optimal connection between the waypoints on the route. This resulted in a modelled 
route network that, while not perfectly fitting the reconstruction, shares important charac- 
teristics with it, especially the presence ofa junction outside the set of known points, as well 
as deviations from the least cost paths that are more in line with the reconstructed route. 
From this we can (cautiously) conclude that the Byzantine road system was constructed 
with the aim of creating optimal connections for the whole micro-regional network. This 
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clearly contrasts with the system of macro-regional connections that characterizes the 
Roman road network as described by the itinerary sources. The results fit with what we al- 
ready know about the Byzantine road system. The new religious“ and military geography of 
the area restructured the micro-regional communication and transport network, using the 
Roman way stations as points of departure. Both the Byzantine settlement pattern and the 
modelled reconstructions suggest that the Roman road itself was not in use anymore. 

As alluded to above, the actual paths used in the Byzantine period may not have been 
optimal with regard to the cost landscape, since their creators would obviously not have had 
access to the full range of cost information, nor to any sophisticated algorithms. An optimal 
topography-based model may be inappropriate and unsuitable for reconstructing past deci- 
sion-making processes (see e.g. Whitley 2002). These, on the one hand, involve a certain 
degree of environmental knowledge (Rockman 2003), and on the other hand concern the 
dynamics of changing land use, settlement pattern and functional or symbolic landmarks 
(military and religious sites; the caravanserai system), which influence the structure of a 
communication system. Depending on the function and aims of the paths and roads, spa- 
tial models can provide multiple optimal paths corresponding to the minimization or maxi- 
mization of specific environmental or cultural criteria, as we did in the case of the visibility 
criterion (see also Howey 2011). However, the choice of which costs to use in the model, 
their values and their respective weights are all subjective, and so results will vary depend- 
ing on these decisions. Furthermore, least cost path modelling will not easily allow us to 
analyse alternative options. The least cost path routines in GIS will give a single optimal sol- 
ution, not the second best or nth-best. Even though k-shortest path algorithms, designed to 
find less optimal solutions, are described in various computer science publications (e.g. 
Yen 1971; Eppstein 1998; Hershberger et al. 2007), they are very difficult to achieve with 
computational efficiency, and have not been implemented in GIS. 

The genetic algorithm described by Frommer and Golden (2007a, 2007b) is in fact 
well suited to address these issues. First of all, the algorithm generates an entire population 
of solutions, not just one. Furthermore, because it has an element of randomness, it may 
not return the same "best" solution each time. Uncertainty in costs and their weightings 
can be reflected in multiple runs of the algorithm. One potentially useful approach would 
be to run the genetic algorithm numerous times with different cost settings and count how 
often a given location is found in the least cost solution (or within a set ofthe top solutions). 
The results of this simulation across various cost value settings would be to assign a prob- 
ability that a path passed through a particular location for each location on the map. 

The study of the development of settlement distribution, density and hierarchy has 
recently been carried out using GIS-based techniques that support the interpretation or re- 
interpretation of the dynamic relationship between supra-regional roads, regional path- 
ways, sites and off-site land use over time (cf. Bell et al. 2002; Fairén Jiménez 2004; Zaksek 


14 For churches as proxies for settlements in the area around Caesarea see Decker (2007, 241). 
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et al. 2008). Other (landscape) theoretical questions still need to be explored, like the in- 
fluence of different environmental costs and sociocultural factors, and landmarks that 
attract movement over the long term (Murrieta Flores 2010), especially where they concern 
the coherence of road networks, their function and the articulation and use by the local 
population of primary (i.e., viae publicae) and secondary networks (regional paths) in the 
Roman period. 

Did the local population of the Cappadocian upland, a territory at the eastern frontier 
of the Roman and then Byzantine Empire, use the structured Roman and Byzantine road 
system? To what extent is the strategic network connecting central places and stations a 
product of the previous regional pattern of sites and infrastructures, and/or related to land 
use and resource exploitation, like mixed farming, agropastoralism and a mining econ- 
omy?» In this respect we also consider it a relevant issue to examine the location and func- 
tion of junctions unifying and subdividing the movement (Gibson 2007) along a diach- 
ronic network. 

The feasibility of exploring these issues is for this case study restricted by the availabil- 
ity of archaeological survey data and fine-grained chronology. However, even for relatively 
well-known regions it is still very much limited by the availability of software solutions. 
Even when using off-the-shelf GIS packages, setting up and executing a multitude of dif- 
ferent least cost path variants and comparing these is a complex task. The current imple- 
mentation of the genetic algorithm for the Euclidean Non-Uniform Steiner Tree Problem is 
suitable for the experimental purpose for which it was applied here, but still needs to be 
substantially improved in speed, flexibility and accessibility to other users before it can be 
used to effect in landscape archaeology. 
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Patricia Murrieta-Flores* 


Developing computational approaches for the study of movement: 
assessing the role of visibility and landscape markers in terrestrial 
navigation during Iberian Late Prehistory 


1 Introduction 


The topic of movement in archaeology has been extensively studied. Research on human 
movement during prehistory has become in archaeology and anthropology one ofthe bases 
for understanding the dynamics of social and economic relationships, technology, social 
identity and territoriality, among many other important themes. Although archaeological 
investigations related to movement have increased in the last decade, the majority have 
usually relied on “static” evidence, that is, on the analysis of the materials or objects that are 
found in specific sites, establishing the relationship between them and their points of ori- 
gin or destination (Branting 2004). In recent years, using spatial technologies, more re- 
search has aimed to investigate movement from a landscape perspective, in which more at- 
tention has been paid to the processes that may have happened on journeys. Some of these 
studies have directly or indirectly analysed the possible factors influencing the decisions 
about which paths to take, the mechanics of movement and the archaeological evidence 
related to it (Llobera 2000; Fairén Jiménez 2004; Cruz Berrocal 2004; Fábrega Alvarez 
2006; Fábrega Alvarez / Parcero Oubiña 2007; Llobera / Slukin 2007; Fiz / Orengo 2008; 
Murrieta-Flores 2010, 2012a; Mlekuz 2010; in the current volume, Lock et al. and Mlekuz 
among others). In the specific case of Iberia, megalithic monuments are among the archae- 
ological elements at a landscape scale that have been linked to potential patterns of move- 
ment, and it has been argued that, besides their symbolic and funerary meanings, they may 
also have been utilized as landscape markers. 

The idea of prehistoric monuments as markers can be traced in archaeology to the 
1970s with the adoption of Saxe's (1970) Hypothesis No. 8 by researchers like Fleming 
(1973), Renfrew (1976) and Chapman (1981). These authors proposed that megalithic 
monuments may have been used as markers with territorial connotations. Nonetheless, 
due to the physical and spatial characteristics of the monuments, other hypotheses have 
also arisen, especially in Iberian archaeology. These theories have pointed out that mega- 
lithic monuments were not only an important element structuring prehistoric landscapes, 
but also that they tend to sustain spatial relationships with historical pathways. This has 
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opened up the question of their possible use as markers, but with a different connotation 
from that of territorial definition. The idea of an association between historical pastoral 
routes and megaliths has been sustained not only by the strong transhumant tradition of 
the Iberian Peninsula, but also by the archaeological record, which points to the practice of 
mixed economies with an important emphasis on mobile pastoralism during late prehis- 
tory (Galán Domingo / Martín Bravo 1991-1992; Ruíz-Gálvez / Galán Domingo 1991; 
Criado Boado 1993; Galán Domingo 1993; Galán Domingo / Ruíz-Gálvez 2001; Murrieta 
Flores et al. 2011a; Murrieta-Flores 2012a; Murrieta-Flores 2013; Wheatley et al. 2010). 
These ideas have resulted in the hypothesis that megalithic monuments may indeed have 
acted as landscape markers, but may also have played a visual role as waypoints in terres- 
trial navigation during the journeys performed by these societies. Although the employ- 
ment of megaliths as landscape markers during prehistory has been implied in many in- 
vestigations, there has never been a specific study on Iberia to test the spatial relationships 
observed in the field. In addition, the concept of landscape markers has been frequently 
used, but never defined in the context of archaeological studies. 

The purpose of this chapter is to investigate the possible role of megalithic monu- 
ments in terrestrial navigation through the development of computational methodologies 
using Geographic Information Systems (GIS). In order to accomplish this, we will test the 
hypothesis of their potential use as waypoints along pathways, using as a case study the ar- 
chaeological evidence available of Copper Age (circa 3100-2100 BC) societies of the moun- 
tain region of western Sierra Morena in Andalucia, Spain. 


2 Landscape markers 


The first task in testing the hypothesis of the use of monuments as landscape markers by 
prehistoric societies is to define what is to be understood as landscape markers. For that, it 
is necessary to take into account that there are many different variables and factors (social, 
economic and environmental) that can influence how people navigate through the land- 
scape (Fairén Jiménez 2004; Murrieta-Flores 2010). Social relationships between groups, 
territorial conflicts, ideas regarding certain spaces or regions, geographical knowledge, 
types of terrain, time expended on travel, velocity of movement, access to basic resources 
while travelling and loads carried are of the greatest importance in the study of movement, 
as they define the routes people take. Along with these and other factors, terrestrial navi- 
gation, or the ability to orient oneself within the landscape, is essential in order to make a 
successful journey. That is, to follow a route, in order to get to a desired destination and to 
return to a point of departure, people have to acquire a series of skills. Research in experi- 
mental psychology has established that wayfinding is one of the most important of these 
abilities, and it consists in the process of defining and following a route from one place to 
another (Golledge 1999). The development of this skill is possible thanks to the creation of 
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cognitive maps, which allow the construction of mental representations of the ideas that we 
hold regarding space and places, the environmental features that can be perceived and the 
spatial connections and relationships that they sustain (Murrieta-Flores et al. 2011b). Thus, 
to carry out a journey, the traveller will have to establish points of origin and destination, 
know the places where he/she needs to deviate, be able to recognize spatial markers and 
understand the complete route within a larger frame (Murrieta-Flores 2010). Therefore, the 
recognition of outstanding features in the landscape, or waypoints, plays a fundamental 
role in human movement (for a complete discussion on the theory and definition of land- 
scape markers see Murrieta-Flores 2011b; 2012b). 


3 Megaliths as waypoints in Iberia 


Although it has been generally accepted that the emergence of megalithic tombs seems to 
be related to the first evidence ofa wish by communities to delimit their space, in the case of 
Iberia they have also been associated with crossing places and markers (Galán Domingo / 
Martín Bravo 1991-92; Andrés Ruperéz 1999; Galán Domingo / Ruíz-Galvez 2001, 264; 
Gómez Vila 2005, 407). Therefore, mobility has been considered one of the main variables 
influencing the location of mounds and megaliths, suggesting that they may have been 
used as waypoints (Criado Boado et al, 1994, 39). 

Observations giving support to this theory have been made since the 19th century, when 
Silvelo (1887, in Gómez Vila 2005, 406) affirmed that megaliths were found along Roman 
paths, and that they were normally positioned on both sides of the road. Later, as a conse- 
quence of early research tracing the Antoninan route, Diaz Sanrujo (1906, in Gómez Vila 
2005, 406) stated that the presence of megalithic monuments was a strong indication of the 
existence of ancient paths, arguing that the most favourable crossing places in a watershed 
would be the ones with megalithic monuments nearby. He also declared that this assertion 
was so accurate that new roads were actually being constructed by identifying these places. 
Later, during the early 20th century and after some decades of discussion of these theories, 
Maciñeira (1943 and 1947, in Gómez Vila 2005, 407) documented 87 monuments along 
40 km of both sides of the Camino Real (royal pathway), and found traces of wheel tracks in 
various sections ofthe excavation along the road. In his report, the author not only argued that 
there was a relationship of intervisibility between the monuments and the road, but also con- 
cluded that, in a landscape with little vegetation, these monuments were in fact acting as way- 
points along the route. More recent spatial analyses have also led to the conclusion that mega- 
lithic monuments were associated with pastoral activities due to their spatial linkage with 
transhumance routes (Criado Boado et al. 1990-91; Criado Boado et al. 1994, 39; Galán Dom- 
ingo / Ruíz-Galvez 2001; Murrieta Flores 2007; Wheatley et al. 2010; Murrieta-Flores 2012a). 

This association seems particularly important in terms of terrestrial navigation, as 
transhumant pastoralism constitutes a strategy that involves the movement of cattle, travel- 
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ling long distances between different altitudes to undertake seasonal exploitation of pas- 
tures in the mountains during the summer and the pastures of the lowlands during the 
winter (Ehlers / Kreutzman 2000, 16). This activity is recorded in Iberia in medieval times 
and, although it has been argued that its practice during prehistory is improbable because it 
was a specific phenomenon produced by the historical circumstances of the Middle Ages 
(Chapman 1979; Walker 1983), it can be suggested that short-distance transhumance (also 
known as trasterminancia) was possible. At the same time, other studies have emphasized 
that the monuments seem not to be related exclusively to historical routes, but also to tran- 
sit zones in general. Research carried out in Aranga, Monfero, La Coruña and Lugo (Gali- 
cia) has established a linkage between mounds, paths and privileged communication zones 
(Criado Boado 1991; Criado Boado et al. 1990-1991, 33-34; Eguileta Franco 1994; Villoch 
Vázquez 2001; Gómez Vila 2005). The region of Cantabria has also been prolific in terms 
of observed associations between menhirs and passageways in watersheds such as Collado 
de los Sejós (Santander) (Bueno Ramírez et al. 1985). A similar association has been ob- 
served for megalithic monuments in the Tajo River Basin (Galán Domingo / Martín Bravo 
1991-1992), where their location seem to be related to vados or fords along the river. In this 
sense, although some of the evidence collected up to now indicates that megalithic monu- 
ments were located at the sides of possible transhumance routes, in other cases it seems 
clear that megaliths were also located at crossing points, perhaps acting as “nodes” or way- 
points along routes, marking or commemorating a specific locus of transit or spots where 
rivers can be crossed. 

The location of the monuments is only one of the qualities observed that is deeply re- 
lated to the theories of why and how megaliths played the role of landscape markers. Monu- 
mentality is another of the characteristics that have led archaeologists to conclude that such 
a use is possible. Prehistoric societies seem to have maximized this effect, constructing the 
monuments in visible or prominent points of the landscape such as hills or crests. In some 
examples it seems that they even enhanced their visibility, using white quartz for their con- 
struction (Bueno Ramírez 1987). However, monumentality as a characteristic of megaliths 
in their use as landscape markers is difficult to assess, as there is a difference or variability 
in the investment devoted to the construction of these monuments in different parts of 
Europe. In terms of movement, differences in monument size probably had repercussions 
in the way a traveller perceived or understood the landscape and, therefore, the role of the 
monuments in mobility. 

In any case, it is beyond reasonable doubt that megalithic monuments were con- 
structed to be noticed by people and to remain as part of the memory of groups and the 
landscape for generations. For these societies, megalithic monuments (funerary or ritual) 
constituted established milestones. As such, not only were they present in the social im- 
agination through their communal and symbolic connotations, but, because of their monu- 
mentality, they also formed part of society's physical reality as principal components of 
their established landscape. In terms of research, although the concept of the use of mega- 
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Fig. 1 | Situation of Sierra Morena in Iberia; western Sierra Morena and all the Copper Age sites recorded. 


liths as landmarks or waypoints has been explored mainly from a theoretical point of view, 
it remains necessary to implement robust methodologies in order to establish whether dif- 
ferent dynamics in their use can actually be observed. As a case study, we selected one of the 
most important mountain regions in Iberia. 


4 The case of western Sierra Morena 


Sierra Morena is a mountain range which constitutes the zone of contact between the Span- 
ish Central Plateau and the Betic Depression and which crosses the northern extreme of 
Andalusia, covering 400 km from east to west (Figure 1). It is located between two of the 
most fertile agricultural river valleys of southern Iberia (Gualdalquivir and Guadiana), and 
it has been characterized by its richness in terms of mineral resources, but also by its ex- 
treme poverty in terms of agricultural potential. In recent years, this region has been the 
subject of a range of diverse projects aiming to understand the settlement patterns and the 
megalithic phenomenon (Pellicer Catalán / Hurtado Pérez 1980; Piñón Varela 1987; 1988; 
Martínez Rodríguez / Lorenzo Gómez 1992; Pérez Macías 1983; 1986; 1987; 1994; 1996; 
1996b; 1996c; 1999; 2010; García Sanjuán et al. 2002; García Sanjuán / Vargas Durán 


DEVELOPING COMPUTATIONAL APPROACHES FOR THE STUDY OF MOVEMENT 


103 


2002; García Sanjuán et al. 2004; García Sanjuán et al. 2006; Hurtado Pérez et al. 2011). 
Thanks to these studies, a more complete picture of the societies that inhabited this moun- 
tain range has been acquired. 

The region of western Sierra Morena (comprising Sierra de Huelva and Sierra Norte 
de Sevilla), has come to be of particular interest due to the fact that societies during Late 
Prehistory, and especially during the Copper Age (3100-2100 BC), seem to have had not 
only a different kind of access to resources compared with those in the valleys, but also an 
alternative social organization. According to the studies of settlement pattern, during this 
period there was a general growth of population, and new settlements were also established 
in places where there were none before and different kinds of habitats were apparently 
used, ranging from larger settlements to apparent seasonal camps. On a regional scale, 
western Sierra Morena acts in many ways as a kind of border area with respect to its sur- 
roundings. In terms of its material culture, and mainly in its pottery, the Chalcolithic sites 
of western Sierra Morena seem to share many characteristics with those in the Guadalqui- 
vir Valley, eastern Andalusia and the southwest in general (Pérez Macías 1994, 4). However, 
its difference from the neighbouring regions seems to be more related to the particular en- 
vironment of the sierra; this is reflected in the subsistence activities, which are closely tied 
to mobile herding, due to the almost non-existent potential for agricultural activities, in 
contrast to the mainly agricultural character of the adjacent valleys. 

From the 388 sites recorded by various projects in the area of the western Sierra 
Morena, 113 megaliths and megalithic cists have been included in the study. Only a few 
menhirs (4) and sites with rock art (2) have been recorded; although it is known that they 
are not the only ones in existence, they are the ones with exact co-ordinates available. 


5 Assessing the use of megalithic monuments as markers through visibility analyses 


As explained before, several theories to explain the observed empirical association between 
megalithic monuments and pathways have been discussed in recent years. In the specific 
case of western Sierra Morena, pastoralism probably started in the Neolithic, and the need 
for pasturage may have been one of the main motors for movement. Recent research has 
concluded that late prehistoric societies of this region probably practised short-distance 
transhumance and that their symbolic sites were deeply related to their dynamics of move- 
ment (Murrieta-Flores 2007; 2012a). These studies have also recorded a strong spatial re- 
lationship between monuments and herding pathways, pointing to the possible role of these 
sites as stations in pastoral orbits. As important spatial markers, one possible explanation 
for this observed relationship is that monuments could have served as waypoints for terres- 
trial navigation during the transhumant cycles. As explained before, megalithic monu- 
ments offer all the characteristics desired in a landscape marker, not only through their dis- 
tinctive shape and size, but also because of their deep symbolic and social importance. 
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To test if they could have been used as waypoints we can make use of GIS visibility 
analysis. As observed above, there are many factors that may influence how humans move. 
However, one of the main and most important prerequisites for navigating while travelling 
is the use of vision. Visibility becomes a critical factor in terrestrial navigation because it is 
through vision that we can identify the markers or waypoints (including celestial bodies) 
that will allow us to take decisions regarding the next direction to take while travelling. In 
this sense, in order to investigate if the megalithic monuments could have acted as way- 
points, we need to establish if they were actually visible along possible travelling routes 
rather than merely close to them. The visibility of the monuments becomes essential as, in 
this case, the possible spatial relationships between waypoints and travellers (and therefore 
the cognitive maps that are constructed) are generated in a dynamic way as the observer 
moves along a route. The specific aim of the experiment presented in this chapter is thus to 
investigate a possible visual structuring of the landscape in relation to the monuments, by 
means of a visibility study that analyses the relationships between natural corridors and the 
monuments as the traveller moves along the paths. 

Although this seems straightforward in a sense, it is necessary to take into account two 
specific considerations in our case study. Firstly, due to the observed close spatial relation 
that the monuments sustain with natural corridors (Murrieta-Flores 20124), the visibility 
patterns that we may detect could be product of a cross-correlation. That is to say, they 
might be visible from the corridors because they are close to them, or, on the other hand, 
they might be close because from those positions they are visible. Both connotations have 
important implications for the interpretation of the patterns left by these societies. In the 
first case, it could imply that the monuments were built along the paths for practical rea- 
sons, like easy transport of building materials, for instance; in the second case, it could 
imply that the monuments were located near to the pathways with the purpose of being vis- 
ible, due to symbolic or traditional activities like ceremonial processions or marking 
stations within herding orbits. Although this is a complex question that will not be solved in 
this study, the following analyses provide some pointers towards fresh ideas. The second 
consideration that needs to be taken into account is that many aspects of terrestrial navi- 
gation are deeply related to movement. While it is possible to confirm through a GIS-based 
visibility analysis whether the monuments are visible from the routes or not, mapping the 
visual changes that occur within the landscape in those terms will not necessarily reveal 
specific perceptions of the landscape while travelling.’ Nevertheless, features created pur- 
posefully for terrestrial navigation usually tend to be strategically located to be as visible as 
possible. In this sense, we can not only explore how constant the view of a monument is 
within the rate of visual change along a route, but we can also identify those locations on 
which one or more features have a visual impact. This can help us to understand the visual 


I For an extended discussion of this and other difficulties such as the role of vegetation, etc., see Murrieta- 
Flores, 2012b. 
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Fig. 2 | Moran's I-Z Score of the megalithic monuments. 


structure of the megaliths from corridors, allowing us to detect whether there was an inten- 
tionality behind the location of these monuments that could be related to both visibility and 
movement. 


6 Constructing visibility methodologies: 
Parameters and technical development of the analyses 


As said above, visibility is subject to important factors related to scale, such as distance. Be- 
cause visibility decays with it and due to the fact that to navigate people have to be within a 
“close” range of a waypoint in order to see it, a sub-regional scale was chosen to carry out 
this experiment. Following the recommendations made by Wheatley and Gillings (2002, 
209) regarding the selection of study areas, it was decided to perform spatial statistics to di- 
vide the entire region into groups of monuments. 

In order to select the study sub-regions, a Kernel Density Analysis was carried out. 
This method is used to calculate the magnitude per unit area of points, identifying possible 
clusters. In this analysis it was considered of importance to establish an appropriated 
neighbouring search radius, because it influences how the kernel density is performed. In 
order to do that a spatial autocorrelation test (Moran's I) was carried out, establishing a stat- 
istically significant radius to carry out the kernel density. From the results obtained in the 
Moran's I test, the peak of the Z score (which measures how intense the cluster is) was 
chosen and used as search radius for the Kernel Density analysis (Figure 2). The output in- 
dicated five major clusters. In reality the Aroche area (western cluster) shows a more con- 
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Fig. 3 | Result of the Kernel Density for the definition of clusters. 


sistent concentration in its northern part, and for this reason, as well as for analytical pur- 
poses, it was decided to derive a sixth group (Figure 3). 

In this manner, the sub-regions established are: Almadén de la Plata (Group 1), Castillo 
de las Guardas (Group 2), Aracena (Group 3), Zalamea la Real (Group 4), Aroche (Group 5) 
and Santa Barbara de Casa (Group 6). It also has to be noted that the region of the Murtigas 
River could not be considered within this analysis due to the fact that there is only one 
monument recorded in that area (Figure 4). 

For the visibility analyses and in order to approximate more accurately human vision, 
the vertical angle range considered was between +10 and —15 degrees, as suggested by the 
results of some physiological studies consulted (Ogburn 2006). In terms of the visibility 
range, while in optimal conditions it has been established that a range between 23-39 km is 
acceptable (Ogburn 2006), in management and archaeological studies the visual ranges 
used have been more conservative (15-18 km and 3-6 km respectively) (Wheatley 1995; 
García Sanjuán et al. 2006). We decided to use a threshold of 3 km, assuming a theoretical 
maximum viewshed of 28.2 km?. 

As said above, the idea behind this experiment was to test if the megalithic monu- 
ments could act as waypoints for terrestrial navigation. According to the theory developed 
on landscape markers (Murrieta-Flores et al. 2011b; Murrieta-Flores 2012b), to fulfil this 
condition, monuments are expected to present one or more of the following character- 
istics: 

(1) they should tend to be considerably visible from the natural corridors or passageways, 
(2) the dominant direction of their fields of view should be those orientated towards the 
natural corridors, and 
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Fig. 4 | The 6 groups defined for the visibility analyses. 


(3) they should be located in a distance range from the corridors that makes them more 
likely to be visible and therefore used as waypoints. 

In order to test these ideas, three different spatial analyses were developed. 

In the case of the first characteristic — that the monuments should tend to be highly visible 
from the natural corridors or passageways — a binary viewshed analysis could be used to 
simply test whether monuments are visible or not from natural corridors.? However, if they 
are also visible from many other locations within the landscape, their visibility may not be 
regarded as a determinant attribute in their role as waypoints for terrestrial navigation. This 
is not to say that in order to function as landscape markers, the monuments should only 
have to be visible from the routes. However, if they acted as waypoints, a closer visual rela- 
tionship between monuments and routes could be expected than between any random 
point in the landscape and routes. In order to test this, we conducted viewshed analyses 
and documented the parts of the natural corridors from which the monuments are visible 


2 These corridors were established in a previous study that also used GIS methods; the full documentation of 
their calculation can be consulted in Murrieta-Flores 2012. 
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Fig. 5 | Segments of the natural corridors from which a particular monument is visible (in red). 
Dolmen de Gabino in Almadén de la Plata. 


(Figure 5). The same was done with a sample of random points, and later both cases were 
compared through a statistical test of significance. The purpose of this was to establish if 
there is a difference between the two samples. The idea behind this is to investigate 
whether the monuments occur in locations that are more visible from the natural corridors 
than would be expected by chance alone. 

The test chosen was the Kolmogorov-Smirnov test and it involves the comparison of 
the cumulative frequency distributions of the two independent samples. In this case the null 
hypothesis set for all groups was: Ho= The megalithic monuments are randomly distributed 
with respect to their visibility from natural corridors. Following the method described, a total 
of six tests were conducted at a significance level of 0.05, assessing each of the sub-regions 
established. 

These procedures were set as the base analysis in order to explore whether we can ob- 
serve significant relationships between the visibility of the monuments and the corridors. 
In addition to this, with the record of the segments of corridors from which the megaliths 
are visible, it was possible to explore certain characteristics between these parts of the land- 
scape and the monuments. In this manner, profiles of the visual exposure to monuments 
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from corridors were also produced from the visibility analysis. This was aimed at graphi- 
cally visualizing the continuity in the visibility of monuments along the corridors, and to 
quantify in each of the sections how many megaliths might have been visible at the same 
time (Figure 6). 

In the case of the second characteristic to be explored - that the dominant direction of 
the monuments' fields of view should be those orientated towards the natural corridors — 
an analysis investigating the possible dominant directionality in the fields of view of the 
monuments was carried out to assess whether these are in places that favour or promote 
their visibility from the routes or, in other words, whether there is a pronounced directional- 
ity in their fields of vision that promotes a visual relationship with natural corridors. The 
method followed here is the one proposed by Wheatley and Gillings (2000). The analysis 
consists in the calculation of a simple binary viewshed of the viewpoint (the monuments in 
this case), the calculation of its Euclidean Direction, then the reclassification of this raster 
into directional zones (N, NE, E, SE, S, SW, W and NW). Finally, the binary viewshed is over- 
laid with the reclassified Euclidean direction in order to quantify the dominant directional- 
ity in the viewshed (Figure 7). 

A table with the results of this analysis was created for each group in order to quantify 
the predominant directions of view of all monuments. Once the dominant directions for 
each megalith were identified, it was examined whether a natural corridor passed through 
it or not. This was used to calculate the percentage of monuments whose main direction 
of view is oriented towards a natural corridor. In order to investigate if the coincidence 
observed between them could be a matter of chance, a 1-sample Chi-square statistical test 
was also conducted. The purpose of this was to identify in a robust way any possible pattern 
of visual structuring that might indicate a preference for locating the monuments in places 
that favoured visibility towards and/or from corridors. In this case, a null hypothesis was 
formulated as: Ho- The natural corridors are equally distributed across dominant and non- 
dominant fields of view of the monuments. The test was carried out at a significance level of 
0.05. 

To test the third characteristic mentioned — that the monuments be located in a dis- 
tance range from the corridors that makes them more likely to be visible and therefore used 
as waypoints — a Higuchi viewshed analysis was conducted. However in this case, instead of 
using the distance index proposed by Higuchi (1988), we created particular ranges for our 
study areas that take into account field observations and the experience of Llobera’s study 
(2007b). Higuchi’s example was adopted for this experiment in terms of the methodology 
followed to establish distance ranges (measuring the angles that a particular object oc- 
cupies in a field of view, and deriving the distances from this). However, instead of using as 
reference a natural element such as trees, we decided to take an archaeological element, in 
this case the megalithic monuments. This was thought to be more appropriate, because we 
are interested particularly in the effect of distance on the visibility of monuments and not 
on the general landscape. 
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Fig.6 | Example of the profile of visual exposure to monuments from a route in Almadén de la Plata. 
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Fig. 7 | a) Simple binary viewshed; b) Calculation of Euclidean Direction; c) Reclassification of Euclidean 
Direction; d) Overlap of viewshed and direction; e) Graphic showing the predominant direction of the viewshed. 
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Fig. 8 | Procedure to generate a Higuchi viewshed: a) Binary viewshed; b) Calculation of Euclidean Distance; 
c) Reclassification to the index established; d) Combination of both. 


The method proposed by Wheatley and Gillings (2000) was followed, but particular dis- 
tance ranges were adopted in this case. The analysis consisted in the calculation of the 
viewshed from the point of interest and its classification into meaningful distance ranges. 
The classification used in this case is based on a field experiment conducted by looking at a 
particular monument and its integration into the general landscape. This was done in a 
similar way as Higuchi did in the case of trees, and Llobera (2007b) in the case of monu- 
ments. Taking a monument with a height of 2.5 m, the visual angle that occupies the monu- 
ment in the field of view of an observer was measured, establishing alpha (a). The distance 
ranges were calculated according to the field observations, taking into account that the 
closer the observer, the wider the angle (Figure 8). The method followed has been described 
elsewhere (Murrieta-Flores 2012b) and it is sufficient to say that the angles were measured 
with an inclinometer. In this manner, the ranges were defined as follows: 

Short Range: a larger than 0.3°. Activities closely related to the monument probably 
happened at this range. At a larger scale the monuments act as reference points in the land- 
scape, they are very well defined and their shape can be easily identified. A person can be 
recognized effortlessly in their surroundings. 
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Fig. 9 | Comparison between the lengths (m) of the path segments from which the monuments 
or random points are visible. 


Middle Range: a from 0.08? to 0.3". The monument can still be a reference point but is 
likely to appear more integrated into the wider landscape. 

Large Range: a less than 0.08. Although the monument may still be visible, it may 
start to fade into the surroundings. How easily it is recognized will probably depend on the 
actual shape of the monument or its contrast with the surrounding landscape. 

In this sense, the maximum distance in the Short Range will be 478.01 m, in the case 
ofthe Middle Range it will be 1,798.56 m and the limit ofthe Large Range will be defined by 
the distance established in the actual visibility analysis, which in this case is 3,000 m. 


7 Results 


The results obtained from the statistical tests were very interesting in terms of observable 
patterns. From the six tests (one per group) conducted to investigate if megalithic monu- 
ments are randomly distributed with respect to their visibility from natural corridors, in 
three cases the null hypothesis was rejected at a significance level of 0.05. While the 
groups in this situation were Almadén de la Plata (Group 1) (25 monuments), Castillo de 
las Guardas (Group 2) (22 monuments) and Aroche (Group 5) (18 monuments), whereas 
for Aracena (Group 3) (16 monuments), Zalamea la Real (Group 4) (29 monuments) and 
Santa Bárbara de Casa (Group 6) (9 monuments), the null hypothesis could not be re- 
jected. 
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Although this appears at first sight to be an inconclusive result, delving into the results 
of the analysis we see that in all cases the megalithic monuments are on average more vis- 
ible from the natural corridors than from random points in the landscape. This is demon- 
strated by comparing the visibility of the random sample and the monuments, where it 
can be observed that the length of segments from which they are visible is more extensive 
(Figure 9). 

Even in the groups where the null hypothesis could not be rejected, the megalithic 
monuments were visible for longer sections than the random points. In the case of the 
Aracena group, while the total length of segments from which the random points are visible 
is 13,951.293 m, the total length from which the megalithic monuments are visible is 
70,550.484 m. In the case of Zalamea la Real, the same phenomenon is observed, where 
the total length of the corridors from which the random points are visible is 6,040.666 m, 
and 13,303.114 in the case of the megalithic monuments. For Santa Bárbara de Casa the dif- 
ference is also large, with a total length of 2,120.460 m for the visibility of random points, 
and 13,184.366 m for the monuments (Figure 10). This means that in general there is a 
larger visual exposure of the megalithic monuments along the natural corridors than other 
points in the landscape, even if we cannot say that the relationship is statistically significant 
at a 0.05 level. 

This is obviously connected to the already proved spatial relationship between monu- 
ments and corridors and it could be thought of as an indication of their use as waypoints. 
However, if the megalithic monuments were acting as such for travellers who were not 
necessarily familiar with the landscape, and their locations were intended as terrestrial 
guidance, they should have acted as a relatively regular system. In other words, along the 
route at least one of the monuments would be visible at all, or almost all, times. This is 
related to the way terrestrial navigation is accomplished, where in order to follow a route we 
need to be able to connect spatially the diverse markers on the route. In order to test this, 
the visual exposure to the megalithic monuments along the routes was graphed, measuring 
for how long segments have a visual exposure to the monuments in each group. For this, 
the number of monuments that are visible each 10 m along a path was recorded, generating 
a visual graphic of its exposure to them (Figures 11 and 12). 

Although in this chapter only some examples are presented, it can be noted that, in 
most of the groups, visual exposure to the monuments along some of the routes was de- 
tected. Nonetheless, there tend to be considerable *gaps" in the visual exposure from 
routes. This contradicts the notion of a visual network. However, a positive outcome is that 
the analysis allowed the identification of sections on the pathways that seem to have 
prompted the visibility of a large number of monuments. This can be regarded as signifi- 
cantin terms ofthe visual impact that megaliths and other prehistoric monuments (such as 
Bronze Age stelae) may have had within the landscape (García Sanjuán et al. 2006; García 
Sanjuán 2011), and the establishment of symbolic places possibly related to specific groups. 
In this manner, the continuous visibility of these monuments from the paths emphasizes 
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—— MONUMENT NOT VISIBLE 


Fig. 10 | Segments of the natural corridors from which one or more megalithic monuments are visible: 
a) Almadén de la Plata, b) Castillo de las Guardas, c) Aracena, d) Zalamea la Real, e) Aroche and f) Santa Bárbara 
de Casa. 
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Visual Exposure 


Fig. 11 | Cumulative visual exposure to megalithic monuments from natural corridors of Almadén de la Plata and 
the longest gaps in visual exposure between (a) the San Bartolomé and Palacio groups and (b) El Chaparral III and 
Dolmen del Esparragal. 
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Visual Exposure 


Fig. 12 | a) Segments in the natural corridors from which the megaliths are visible. b) Cumulative visual exposure 
to all megalithic monuments from natural corridors of Aroche. 
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Fig. 13 | Visual exposure of the megalithic groups of San Bartolomé and Dolmen de Palacio III from a segment 


of the natural corridor that runs along the transhumance route of Cordel del Pedroso. 


their significance in terms of the creation of a funerary and symbolic landscape, where the 
visual accessibility of the monuments seem to have had great importance. This is sup- 
ported by the results obtained by the statistical test already mentioned, which was con- 
ducted in order to investigate directional viewsheds. This analysis showed a significant re- 
lationship in most of the groups and that the majority of the dominant directions in the 
fields of view from monuments are related to a natural corridor or to a transhumance route 
(Figure 13). However, in terms of terrestrial navigation this could be interpreted in a differ- 
ent way. Although the megaliths in many cases are quite well distributed along the main 
corridors, many ofthem are clustered (Figure 13). Whether or not they are all contemporary, 
in neither case would it have been necessary to have several waypoints very close to each 
other if one of them was already visible from the corridor. In this manner, the adjacent lo- 
cation of the monuments is perhaps more related to a symbolic recognition of particular 
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Fig. 14 | Percentages per group of monuments that present a field of view with a dominant direction towards 
natural corridors. 
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places that acted as mnemonic markers for a number of groups, families or households, 
but not necessarily as waypoints. 

Regarding the directional viewsheds, the general results were very interesting. Al- 
though in the first statistical analysis the relationship between general visibilities and the 
corridors in some ofthe groups could not be confirmed, the tests exploring the specific domi- 
nant directions in the fields of view of the monuments and their possible connection to cor- 
ridors was statistically significant. This was the case in all groups with the exception of Santa 
Bárbara de Casa, where the test could not be carried out due to the small sample. A closer 
look at the individual results of this analysis confirms that in 5 of the 6 groups, more than 
8096 ofthe monuments exhibit fields of view in which the dominant directions are the ones 
that coincide with the natural corridors. Again, this is the case in all groups with the excep- 
tion of Santa Bárbara de Casa, where the dominant directions are related to a natural corridor 
in only 66.6696 of the monuments. While in Almadén de la Plata 80% of the megaliths 
present a relationship between predominant direction of view and corridors, for Castillo de 
las Guardas the figure is 86.36 %. In the case of Aracena the coincidence is 81.25 96, for Za- 
lamea la Real 80.9596 and for Aroche 94.44% (Figure 14). These results point to the idea 
that, towards the Copper Age, the majority ofthe monuments within the sierra were built in 
locations from which natural corridors were visible. In fact, the results suggest that locations 
with a dominant view oriented to the corridors were preferred (Figure 15). 

The modified Higuchi viewsheds reaffirmed the observations made previously, which 
pointed to a close visual relationship between monuments and corridors. It was of great in- 
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Fig. 15 | Notice the coincidence between the dominant view direction of La Portilla (above) and Pasada del Abad in 
Aroche (below) with natural corridors. 


terest to observe that the majority of the natural corridors fall in what we defined as a visual 
Middle Range (Figure 16). This seems to be consistent also with the visibility of historical 
routes, where, as in the previous case, the majority also fall in this range. Nevertheless, 
there seems to be a closer relationship between monuments and natural corridors than be- 
tween monuments and historical routes (Figure 17). It is also worth mentioning that, al- 
though the number of natural corridors that fall in the visual Short Range is lower, in the 
majority of the groups there is still a very high coincidence between corridors and monu- 
ments where in up to 66.6696 of the cases a natural corridor passes through the shortest 
range. This is relevant because a close visual proximity to the corridors implies a greater 
chance of monuments being visible, and therefore that they could be used as waypoints. 
However, a more consistent relation to the Middle and Long Ranges is observed, and this 
seems to suggest that there was an intentional relationship of proximity between them. 
This is also supported by the fact that even in flatter and broader landscapes like the Viar 
Valley (Almadén de la Plata, Group 1) in which the options about where to walk are very 
wide, late prehistoric societies of this region chose to locate their monuments at close visual 
ranges to natural corridors. As said, this is observed in the majority of the cases and this re- 
lationship could be associated with the way herding trips were conducted and the way that, 
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Fig. 16 | Results of the Higuchi Viewsheds: a) Percentage of coincidence between natural corridors and visual 
distance ranges established and b) between historical transhumance routes and visual distance ranges. 


Fig. 17 | Notice the coincidence between the dominant view direction of La Portilla (above) and Pasada del Abad in 
Aroche (below) with natural corridors. 
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as places to visit during these journeys, the visibility of symbolic places in relation to corri- 
dors became important. It is very interesting to note that other research (Galán Domingo / 
Ruíz-Galvez 2001; Murrieta-Flores 2012a; 2013) has suggested that megaliths may have 
also served as landmarks at fords and mountain passes, and this observation seems to be 
strengthened by this analysis. 


8 Conclusion 


The assessment of the use of megaliths as waypoints has proved a difficult task as the mod- 
els are relatively simple and the variables to take into account are complex. However, it can 
be said that the experiments were successful in the identification of a series of interesting 
patterns. The relationships observed between dominant directions in the fields of view and 
corridors in the majority of the groups, and the high coincidence between general visibil- 
ities and pathways, point to the idea that monuments were effective landscape markers. 
However, they seem not to have been built with a focus on terrestrial navigation. Itis certain 
that there is a visual relationship between the monuments and natural corridors, and 
people familiar with this landscape may have been able to locate them easily. Although 
there is no single corridor where monuments could be regarded as part of a “complete” vis- 
ual system in terms of continuous visibility, it must be also taken into account that many 
other features in the landscape may have acted as intermediate waypoints, such as charac- 
teristic topographical features or trees. In this manner, it is acknowledged that, with in- 
formation such as verbal accounts, foreigners too may on some occasions have been able to 
use them as points of reference while navigating through the landscape. In those terms, 
they may have been used as waypoints too. 

What seems to be certain is that the monuments marked vistas for people moving 
along the corridors. This seems to be confirmed by the fact that the megaliths are seen 
more often from the corridors than are other points in the landscape, and because the 
lengths of the segments from which they are visible within the paths are considerably 
larger. This outcome is reinforced by the results from the directional viewsheds. The fact 
that the dominant directions in the fields of view of the megaliths are the ones oriented to- 
wards the natural corridors suggests that the visibility of the monuments from them and 
the other way around was regarded as important in some sense. In this manner, the reason 
why these two features were visually and spatially connected could lie in the symbolic role 
that they probably played within pastoral orbits (Murrieta-Flores 2012a; 2013). While carry- 
ing out these cyclical movements, mobile herding societies visit places they regard as sig- 
nificant in terms of symbolism and memory. These places are normally related to their sea- 
sonal cycles, and therefore they are located within their areas of movement. Some of them 
are visited through several generations and different groups may also have places in com- 
mon where they can meet to exchange news, celebrate festivities, conduct commerce or 
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marriages or simply spend time with each other. Furthermore, the link between monu- 
ments and movement may be connected to the idea of permanence, helping to mark the 
pastoral orbits of different groups. In this sense, the association of monuments and corri- 
dors may find its explanation in the symbolic meaning of these sites and their use as land- 
scape markers over generations. The possible meaning of the megaliths as the “houses” of 
the ancestors could have played an important role in the visits established during the pas- 
toral circuits, where the long use of sites like Palacio III in Almadén de la Plata or Valdeli- 
nares in Aracena may constitute testimony of this. 

Although it is clear that distinctive locations with particular characteristics were 
chosen to build the megaliths in this region, it is also interesting to observe that in some 
cases they are close also to the sources of material from which the orthostats that compose 
them were taken. This is the case with monuments like Valdelinares in Aracena (Romero 
Bomba 2001, 377). There seems thus to be a component of economy and practicality in the 
choice of these places, but the importance of this factor in determining the selection of the 
sites can only be explored further through detailed study of the sources of building materi- 
als. In addition, it could be said that in many cases the monuments seem to be located near 
rivers and fords. While as part of a herding itinerary these locations would have many con- 
venient qualities in marking sources of water and river crossings, there could also be a sym- 
bolic, where the monument commemorates not only ancestors, but also the importance of 
movement in their way of life. It is of great interest to observe that while settlements seem 
to be seasonal, the only permanent features are precisely these monuments, which, being 
located along natural corridors, could have celebrated the cycles of life and the temporal 
character of their actions, in which movement knit together all meanings and represen- 
tations. 

Therefore, although monuments may not have acted as waypoints in the same manner 
as the Roman miliarium did, it is very possible that they acted as landmarks in the same way 
that other, later societies would recognize the architecture of their own times. Just as a pil- 
grim knows her parish and could identify a Christian church even if far away from home, it 
is possible that prehistoric communities utilized monuments as markers of important 
places within their landscape and, even if they may have occasionally served as a point of 
reference for terrestrial navigation, it is unlikely that all ofthem functioned as such, so as to 
form waypoints within travel "itineraries." 

The location of important symbolic places throughout cultures and times has usually 
been attached to particular events, beliefs or social memories. Although nowadays the spe- 
cific meanings of megaliths escape us, it is interesting to observe that we can still recognize 
in their locations aspects that seem to have been of importance. The constant observation of 
a spatial relationship between megalithic monuments and natural corridors could be some- 
thing specific to the mountain environment of Sierra Morena. However, this connection 
may be related to the dynamics of these societies' daily life. It is possible that hunting and 
herding itineraries were formed over time, where contacts with different groups made 
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possible a “unified” idea in the symbolism of a megalithic landscape, which at the same 
time took a distinctive and particular character in each region. This seems to be the case in 
the western Sierra Morena, where the megalithic phenomenon as well as other archaeologi- 
cal evidence is testimony of the contact between regions. Nevertheless, these contacts seem 
not to be very extensive and the general day-to-day dynamics probably had a more local na- 
ture. This does not mean that they did not actually have extensive contacts with other re- 
gions. Several stone tools and green stones such as variscite have been found in other 
Copper Age sites and they may have their origin in western Sierra Morena (Linares Catela et 
al. 2010; Linares Catela / Odriozola Lloret 2011). 

As said before, the differences observed in settlement patterns, as well as in sizes and 
styles of construction of the megalithic monuments between regions such as the Guadal- 
quivir Valley and Sierra Morena, seem to support the idea not only of a differentiation in 
demographic numbers and access to resources, but possibly also in ideology and the con- 
struction of a symbolic landscape. It even seems possible to differentiate between regions 
within western Sierra Morena. However, what they may have in common is that their mega- 
liths were built by groups formed from several families and the locations selected were al- 
ways near to natural corridors and possibly related to their herding itineraries. In this 
manner, during hunting and herding journeys, the places once regarded as special were 
chosen to become sites with an enduring character within social memory, where the mega- 
liths would effectively become a physically permanent station among different groups, 
marking the space symbolically. 

Later, towards the Bronze Age, the dynamics seem to change with the establishment of 
more permanent settlements like El Trastejón (Aracena). The de-monumentalization of fu- 
nerary practices and the change to individual tombs marked a very significant trans- 
formation of the monuments in terms of social meaning. Nevertheless, in terms of cosmol- 
ogy, these structures and places were still regarded as important in western Sierra Morena. 
This is evident in monuments like Valdelinares (Aracena) which, due to its dimensions and 
construction system, could be regarded as related to a later period during the Copper Age, 
and closer to the megalithic cists. The reuse through many generations of these monu- 
ments and sites, like Valdelinares, Palacio III (Almadén de la Plata), El Pozuelo (Zalamea la 
Real), El Trastejón and La Papúa (Aracena), are testimony of their role as markers and sites 
of power within the landscape. However, although there seems to be continuity in the fu- 
nerary patterns, where monuments were reused or new ones were still located very near to 
natural corridors, this coincidence decreased. Gradually, groups transformed into societies 
where individuals had more weight, a phenomenon that emerged along with an increasing 
social inequality. The apparent establishment of a stronger mobile tradition in terms of 
extended networks, and probably the increase in commerce and exchange, seem to be the 
reasons why we can observe a closer and growing spatial relationship between settlements 
and corridors towards the Early and the Late Bronze Age. With the presence of permanent 
settlements and the development of a society with more social inequality, both settlement 
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and funerary spatial patterns changed. The way of carrying out the pastoral orbit probably 
did not involve families anymore, but one or a few members, and, while the funerary sites 
were possibly part of the itineraries, the settlements seem to have taken on a more import- 
ant role within them (Murrieta-Flores 2013). 
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Abstracts 


D. Mlekuz, Approaching past movement: time to move a little further? 


The paper discusses the ways to approach the past movement using GIS in archaeological 
landscape studies. Through movement we actively participate and create the world around 
us. By moving we bind things, places, persons together in a process of continuous becom- 
ing. Thus movement creates time, space and material world. 

However, in archaeology we are dealing with the static, material residues of this past. 
The archaeological landscape is network of places, which were reproduced as “collapsed 
acts” of persons moving and performing routines in space and time, themselves being 
changes in the process. They are time and movement materialised or 'fossilised' in the 
landscape. Paper presents some new tools and approaches for understanding past land- 
scapes, especially tools for the analyses of proximity, mobility, accessibility and interaction. 
GIS procedures were developed to help us understand the limits time and space imposes in 
interaction between people and help us reveal deeper logic of social processes in landscape. 
Approaching past mobility also forces us to explicitly address the notion of time in GIS and 
archaeological landscape studies. 


Gary Lock, Mariza Kormann Rodrigues and John Pouncett (Oxford), Visibility and 
movement: towards a GIS-based integrated approach 


It is argued here that movement can only be fully experienced and, therefore, realistically 
modelled through integration with visibility. In this paper we explore a range of different 
approaches and techniques that investigate the characteristics of movement and visibility, 
the relationship between them and, importantly, how these are affected by changing scale. 
Using the ancient Ridgeway track in central England as a case-study, we first explore the 
idea of ‘intentionality’ within movement. This uses different sized windows of analysis to 
model looking ahead' and the idea of moving using waypoints as markers. In a related ap- 
proach we develop Llobera's idea of landscape ‘prominence’ and introduce a new method, 
Fuzzy Cumulative Visibility Analysis (FCVA), which simulates a 'virtual walk' in the land- 
scape by moving from cell to cell and calculating visibility in all directions. Some cells are 
seen more often than others in a continuously changing way and these changes can be ex- 
pressed by a visibility index for each cell resulting in visibility gradients due to cumulative 
effects. This approach is novel in being applicable at a range of scales and is demonstrated 
through a case-study of Neolithic long barrows in central England. The approach to each 
barrow can be modelled in terms of visibility, distance and position to show interesting dif- 
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ferences between basing these analyses on topography, prominence and FCVA. Building 
on this we continue with a re-assessment of the idea of Viewpaths' as first suggested by Lee 
and Stucky in 1998. This uses the notion that visibility is an actual 'cost' in influencing 
movement and is modelled by using forms and combinations of visibility indeces as cost 
surfaces for generating least cost paths, corridors and surfaces. Overall, this developing 
suite of approaches provides a much richer understanding of perceiving the landscape 
through movement and visibility than just viewsheds or cost paths. 


E. Fovet, K. Zaksek, Path modelling and settlement pattern 


This paper describes the contribution of path modelling to the ancient settlement pattern 
study over the long term. The path modelling methodology is a stimulating tool, which is 
complementary to the hierarchical approaches in the landscape archaeology since it con- 
tributes to the understanding of the spatial relation between archaeological sites. The exist- 
ing methodology was enhanced by enlarging the set of path reconstitution parameters (vis- 
ibility) and by modelling in two scales. The proposed model is based on parameters derived 
merely from the relief because its changes should be insignificant even over a long period 
time. 


P. Verhagen, S. Polla, |. Frommer, Finding Byzantine junctions with Steiner Trees 


The Byzantine road system in Cappadocia (Turkey) is known in reasonable detail. Recon- 
structions of the road system by Hild (1977) show the most probable location of the routes 
that were built and maintained by the imperial administration. Unlike the earlier Roman 
routes, these were not only built to establish fast connections between major cities; rather, 
they were intended to connect all administrative and religious centres in the area. A con- 
spicuous element of this connecting network is the fact that road junctions are found in the 
countryside, without any corresponding traces of settlement. Least cost path models trying 
to connect the centres will inevitably fail to identify these junctions, since LCP modelling 
will only create routes between single sources and targets, and not find the optimal network 
solution to connect three or more destinations. Steiner trees are an appropriate mathemat- 
ical solution to deal with this problem; however, applications in GIS are lacking up to now. 
In this paper we would like to present the preliminary results of a case study showing that it 
is technically possible to use Steiner trees in conjunction with cost surfaces. The modelling 
results seem promising, as the resulting connections are closer to the documented ones 
than those resulting from least cost path modelling. However, there are substantial unre- 
solved issues regarding the implementation of this method in GIS. 
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P. Murrieta Flores, Walking through mountains and sacred monuments: insights to the role 
of natural areas of transit during Late Prehistory in South-Western Iberia 


Through time, societies have not only learned how to utilise resources from the environ- 
ment, but they have also acquired specialised knowledge, adapting to it. Mountain regions 
are often very special environments due to their diversity, strategic significance and margi- 
nality. Mountain contexts are of great importance for people living within and neighbour 
areas due not only to their richness in natural resources and their great biodiversity, but also 
because of their peculiar socio-cultural dynamics, whose development has tended to be at 
the margins located often far away from the centres of political power. This marginality can 
be also extended to movement within these regions, where mountain ranges regularly con- 
stitute mighty obstacles due to their natural configuration, playing a central role in strategy, 
commerce and travelling. Movement in these environments is frequently localised and ac- 
commodated to the nature and steepness ofthe terrain. Travelling through natural passage- 
ways and optimal routes such as mountain passes and internal valleys is regular in these re- 
gions, and although the optimal employment of the landscape not always prime in human 
societies, in certain special conditions such as these, the search for optimality can have 
serious advantages. In the case of western Sierra Morena (Spain), its configuration shaped 
both, the very forms of transit through the mountains and the historical routes of com- 
munication that traverse Andalucia. In this paper we will explore a GIS methodology devel- 
oped specifically to identify particular characteristics of the landscape that are relevant for 
movement such as passageways, crossing points and natural areas of transit. In our par- 
ticular case study, the identification of such areas allowed us: (1) to examine the role that 
natural accessibility might had for late prehistoric societies of Southern Iberia in the lo- 
cation of their habitats and symbolic places; and (2) to determine through the examination 
of historical and modern paths, to what extent zones of optimal accessibility have been used 
through time in the establishment of communication routes. 
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Abbreviations 


a.s.]. above sea-level 

DEM digital elevation model 

DTM digital terrain model 

FCVA Fuzzy Cumulative Visibility Analysis 
FOV field of view 

GIS Geographical Information System(s) 
LCP Least Cost Path(s) 

Mp Roman mile (milia passuum) 

OD ordnance datum 
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